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ABSTRACT 

The  objective  of  the  STAR  project  is  to  test  and  space 
qualify  a  new  continuous  cycle  cryogenic  refrigeration 
system  for  cooling  of  sensors  and  electronics  which  is  based 
upon  the  newly  discovered  thermoacoustic  heat  pumping 
effect.  The  new  refrigerator  has  no  sliding  seals,  a  cycle 
frequency  of  about  300  Hz,  and  uses  acoustic  resonance  to 
enhance  the  overall  power  density  and  efficiency.  This 
thesis  is  concerned  specifically  with  the  design  and  testing 
of  the  electrodynamic  transducer  which  is  responsible  for 
the  electro-acoustic  power  conversion.  A  computer  model  of 
the  driver/resonator  system  is  presented  along  with  the 
techniques  for  measurement  of  the  electrical  and  mechanical 
parameters  used  as  input  for  the  model. 

A  final  driver  design  (including  dimensional  drawings) 
utilizing  a  modified  JBL  2450J  neodymium-iron-boron 
compression  driver  and  associated  leak-tight  electrical 
feed-throughs ,  microphone,  accelerometer ,  pressure  gage, 
pressure  housing,  and  resonator  interface  is  provided. 
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I .   INTRODUCTION 

A.   BACKGROUND 
1 .   History 

Thermoacoustics  can  generally  be  described  as  the 
study  of  the  interaction  between  heat  and  sound.  Scientific 
interest  in  this  interaction  is  not  new.  Lord  Rayleigh  (1878 
and  1945)  discussed  various  qualitative  examples  of  the 
production  of  sound  jy  heat.  In  one  of  these  examples  he 
describes  investigations  by  Sondhauss  (1850)  of  an  effect 
observed  by  glassblowers :  heating  a  bulb  of  glass  at  the 
end  of  a  hollow  rod  occasionally  led  to  the  generation  of 
sound.  Faraday  (1318)  demonstrated  that  this  effect  occurred 
with  gases  other  than  air.  Sondhauss  proved  that  the  vibra- 
tion of  the  glass  itself  did  not  generate  the  sound,  but  he 
offered  no  explanation  as  to  what  did.  In  his  description  of 
Sonhauss '  work,  Rayleigh  stressed  the  importance  of  the 
phase  difference  between  temperature  and  particle  motion. 
Rijke  (1859)  built  an  open  tube  with  a  wire  mesh  inside. 
When  the  mesh  was  heated,  the  tube  produced  a  sound  of 
impressive  intensity.  The  functioning  of  this  oscillator  is 
slightly  different  since  it  requires  steady  gas  flow  for  its 
operation.  A  similar  sound  production  effect  was  observed  by 
Taconis ,  et  al .  (1949)  in  hollow  tubes  immersed  in  liquid 
helium   (4.2c^j,r^_   The   Taconis  oscillations   were   addressed 


quantitatively  by  Yazaki ,  et  al .  (1980).  The  work  described 
above  deals  primarily  with  prime  movers--devices  that  con- 
vert a  temperature  gradient  to  sound  energy.  We  call  this 
the  classical  era  of  thermoacoustics . 

Interest  in  thermoacoustics  was  renewed  when  the 
idea  of  the  use  of  acoustical  devices  as  refrigerators 
appeared.  Gifford  and  Longsworth  (1966)  described  the  pump- 
ing of  heat  along  a  surface  caused  by  a  periodic  change  in 
the  pressure  of  the  adjoining  gas.  Such  a  change  can  be 
produced  by  the  oscillations  of  a  sound  wave.  In  their 
experiments  Merkli  and  Thomann  (1974)  explored  the  heating 
and  cooling  effects  on  the  wall  of  a  gas-filled  resonant 
tube.  They  found  that  heat  was  transported  from  a  region 
near  the  velocity  antinode  (or  maximum  value)  of  the  sound's 
standing  wave  to  regions  near  the  adjacent  pressure  anti- 
nodes.  Thermoacoustic  prime  movers  have  also  received  recent 
attention.  Kempton  (1976)  investigated  the  excess  noise  of 
aeroengines  above  that  predicted  by  theory.  He  determined 
that  the  sound  was  produced  by  unsteady  heat  transfer.  Each 
of  these  latter  three  groups  of  exprimenters  used  some 
theory  for  comparison,  but  it  was  mostly  qualitative.  None 
provided  the  complete  quantitative  theory  that  would  explain 
their  experimental  observations. 

The  theoretical  breakthrough  in  thermoacoustics 
field  was  made  by  Nicklaus  Rott  (1969,  1974,  1975,  and 
1980).   His   theory  combined  basic  principles  from   physics, 


thermodynamics,  and  acoustics  to  quantitatively  describe  the 
effects  found  in  both  types  of  thermoacoustic  devices ;  prime 
movers  and  heat  pumps.  Prime  movers  use  a  temperature  gra- 
dient to  create  sound,  as  in  the  effects  discussed  by  Sond- 
hauss,  Rayleigh,  Taconis,  and  Kempton.  Heat  pumps,  or  refri- 
gerators, use  the  oscillating  pressure  of  a  sound  wave  to 
produce  a  temperature  difference,  as  described  by  Gifford 
and  Longsworth  and  Merkli  and  Thomann.  Rott  described  the 
effect  found  by  Merkli  and  Thomann  as  thermoacoustic  stream- 
ing. It  is  this  effect  that  makes  thermoacoustic  refrigera- 
tion possible. 

Inspired  by  Ceperley's  (1979)  traveling  wave  Ster- 
ling Cycle  heat  engine  and  Rett's  quantitative  theory,  the 
team  of  Wheatley,  Hofler,  Swift,  Migliori ,  and  Garrett 
(1982,  1983a,  1983b,  1985,  and  1986)  developed  a  series  of 
thermoacoustic  experiments  at  Los  Alamos  National  Laboratory 
in  New  Mexico.  They  investigated  the  basic  thermoacoustic 
effects  in  both  prime  movers  and  refrigerators  and  compared 
their  experimental  results  to  Rett's  theory. 
2 .   Thermodynamics 

We'll  digress  from  history  here  to  explain  the  ther- 
modynamic distinction  between  prime  movers  and  refrigera- 
tors. This  discussion  follows  Sears  and  Salinger  (1975) .  A 
prime  mover  receives  heat  at  a  high  temperature,  T„  (the  hot 
reservoir),  does  work  on  its  surroundings,  and  rejects  heat 
at  a  lower  temperature,  Tc  (the  cold  reservoir),  as  shown  in 


Figure  I-l.a.  The  first  law  of  thermodynamics  tells  us  that 
the  heat  flow  (Qm)  from  the  hot  reservoir  must  equal  the 
work  done  on  the  surroundings  (W)  plus  the  heat  flow  (Qc)  to 
the  cold  reservoir: 

Qm  =  W  +  Qc.  (I-l) 

The  second  law  of  thermodynamics  states  that  in  every  pro- 
cess the  entropy  of  the  universe  must  either  remain  constant 
or  increase.  The  entropy  change  for  the  isothermal  heat 
transport  to/from  the  reservoirs  is  defined  as  the  heat  flow 
divded  by  the  reservoir  temperature.  For  a  prime  mover  the 
second  law  can  be  written: 

^^M  i  ASo    or    Qm/Th  i  Qc/Tc.  (1-2) 

The   thermal  efficiency  of  the  prime  mover  is  defined  as  the 

ratio  of  the  work  output  to  the  heat  input: 

n  =  w/Qm  =  (Qm  -  Qc)/Qm  =  1  -  Qc/Q«.        (1-3) 
Comparing   equations   1-2  and  1-3  gives  us  a   limit   on   the 

efficiency  based  on  the  temperatures  of  the  reservoirs: 

^  ^  1  -  Tc/Tm.  (1-4; 

This  limit  is  known  as  the  Carnot  efficiency. 

In  a  refrigerator  the  process  is  basically  reversed 
(see  Figure  I-l.b).  Heat  (Qc)  is  removed  from  the  low  tem- 
perature reservoir,  work  (W)  is  received  from  the  surroun- 
dings, and  heat  (Q„)  ±s  delivered  to  the  higher  temperature 
reservoir.  The  expression  for  the  first  law  is  the  same  as 
for  the  prime  mover,  Equation  I-l.  For  the  refrigerator  it 
means   that  the  heat  flow  into  the  hot  reservoir  must   equal 
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Figure  I-l. 


The  Thermodynamics  of  (a)  Prime  Movers 
and  (b)  Refrigerators. 


the  heat  flow  from  the  cold  reservoir  plus  the  work  done   to 

the  system.  The  second  law  for  this  case  is: 

ASh  '£    ASc   or   Qm/Tm  ^  Qc/Tc.  (1-5) 

The   efficiency  parameter  for  the  refrigerator  is  called  the 

coefficient  of  performance  (COP) .  It  is  defined  as  the  ratio 

of   the  heat  flow  from  the  cold  reservoir  to  the  work   input 

to  the  refrigerator,  or: 

COP  =  Oc/W  =  Qc/(Qm-Qc)  ^    Tc/{Tm-Tc).         (1-6) 
The  limit  is  known  as  the  Carnot  coefficient  of  performance. 

The  work  in  this  thesis  deals  solely  with  xeifrigera- 
tors .  Even  though  prime  movers  are  mentioned  from  time  to 
time.  the  primary  focus  of  the  rest  of  our  discussion  will 
be  refrigerators . 

3 .   Overall  Efficiency 

According  to  Hofler  (1936),  the  early  thermoacoustic 
work  at  Los  Alamos  focused  on  experimental  refrigerating 
engines.  The  performance  of  these  engines  fell  short  of 
expectations,  leading  to  simple  experiments  on  basic  thermo- 
acoustic effects  and  a  proof -of-principle  experiment  on 
thermoacoustic  refrigeration.  Hofler  then  applied  the  Rott 
theory  to  the  experimental  systems  and  solved  the  resulting 
equations  numerically.  For  his  doctoral  dissertation  from 
the  University  of  California.  San  Diego,  Hofler  designed  and 
constructed  a  completely  functional  thermoacoustic  i^efriger- 
ator .  He  also  made  accurate  measurements  of  its  thermody- 
namic efficiency,   and  used  this  efficiency  to  make  compari- 


sons  to  the  Rott  theory.  After  receiving  his  doctorate, 
Hofler  came  to  the  Naval  Postgraduate  School  (NFS)  as  a 
post-doctoral  fellow,  and  brought  his  prototype  refrigerator 
with  him. 

The  purpose  of  this  thesis,  in  conjunction  with  the 
work  of  several  other  students,  is  to  modify  Hofler 's  re- 
frigerator design  in  order  to  improve  its  overall  efficien- 
cy and  make  it  suitable  for  space  cryocooler  applications. 
The  basic  design  of  the  thermoacoustic  refrigerator  is  shown 
schematically  in  Figure  1-2.  The  driver  (A),  which  produces 
the  sound  waves,  is  coupled  via  bolts  to  the  resonator  (D) 
via  a  reducer  cone  (B)  and  bellows  (C).  Inside  the  resonator 
is  a  stack  of  plastic  plates  (E)  and  their  heat  exchangers 
(F)  which  allow  the  heat  to  be  removed  from  the  hot  end  and 
absorbed  by  the  cold  end.  It  is  the  interaction  between  the 
sound  waves  and  the  plastic  plates  that  produces  a  tempera- 
ture difference  across  the  plates  and/or  pumps  heat  from  the 
cold  heat  exchanger  to  the  hot  heat  exchanger.  A  brief 
description  of  efficiency  would  be  useful  here.  Figure  1-3 
shows  a  diagram  of  the  chain  of  energy  transformations  which 
occurs  in  a  thermoacoustic  refrigerator.  There  are  three 
conjugate  variable  pairs  (V  and  I;  P  and  U;  and  iT  and  6c), 
and  two  transformations  (electrical-to-acoustical  and  acous- 
tical-to-thermal) . 

We  start  with  electrical  input  power  to  the  driver, 
which  can  be  calculated  (using  r.m.s.  values)  with: 


Figure  1-2. 


Basic  Thermoacoustic  Refrigerator  Showing 

(A)  Driver,  (B)  Reducer  Cone,  (C)  Bellows 
(D)  Resonator,  (E)  Stack,  and  (F)  Heat 
Exchangers . 
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Chain  of  Energy  Transformations  in  a 
Thermoacoustic  Refrigerator. 


P-i-c=  =  I  V  COS  e,  (1-7) 

where  I  is  the  current  into  the  electroacoustic  driver.  V  is 
the  voltage  across  the  driver  terminals,  0  is  the  phase 
angle  between  I  and  V,  and  cos  e  is  known  as  the  "power 
factor."  An  electrodynamic  driver  converts  the  electrical 
power  to  acoustic  power  (the  first  transformation).  Analo- 
gous  with   the  previous  definition  of  the  electrical   input 

power,  the  acoustic  power  is  given  by: 

PJ  =  P  U  cos  ^,  (1-8) 

where  P  is  the  acoustic  pressure,  U  is  the  volumetric  velo- 
city which  equals  the  product  of  particle  velocity  and 
resonator  cross-sectional  area,  and  i^  is  the  phase  angle 
between  them.  The  efficiency  of  this  electroacoustic  trans- 
formation is  given  by: 

ne«  =  t^/Pe^ec.  (1-9) 

In  the  second  transformation,  the  refrigerator  con- 
verts the  acoustic  power  to  a  temperature  gradient  (vT)  and 
heat  flow  iCi)  .  As  discussed  previously,  the  efficiency  for 
this  transition  is  given  as  the  coefficient  of  performance: 

COP  =  (^c/^  .  (I-IO) 

The  overall  thermoacoustic  efficiency  is  therefore  the  pro- 
duct of  the  efficiencies  for  the  two  transformations: 

n-r«  =  T1e«  ♦  COP  =  (^c/Pelec.  (I-ll) 

Improvements  in  the  overall,  refrigerator  efficiency 
can  thus  be  divided  into  two  distinct,  though  interrelated, 
efforts.   The  first  is  the  subject  of  this  thesis:  modifica- 
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tions  to  the  driver  to  optimize  the  electroacoustic 
efficiency.  The  second  effort  involves  the  optimization  of 
the  coefficient  of  performance  and  is  addressed  by  LT 
Michael  P.  Susalla,  USN ,  in  his  master's  thesis. 

B.   SPACE  THERMOACOUSTIC  REFRIGERATOR  (STAR) 
1 .   Motivation 

In  addition  to  improving  overall  efficiency,  our 
experimental  refrigerator  will  be  designed  and  built  with 
the  intention  of  flying  it  on  the  Space  Shuttle  as  part  of 
the  National  Aeronautics  and  Space  Administration's  (NASA) 
Get  Away  Special  (GAS)  program.  As  of  3  February  1988,  there 
is  a  Memorandum  of  Agreement  between  the  Naval  Postgraduate 
School  (NPS)  and  the  Air  Force  which  provides  funding  for 
the  space  flight  and  which  assigns  the  NASA  payload  number 
G-337  to  this  project.  The  reason  for  testing  the  thermo- 
acoustic  refrigerator  in  space  is  that  the  immediate  appli- 
cations of  the  STAR  are  expected  to  be  the  cooling  of  elec- 
tronics, high  Tc-  superconductors,  and  infrared  detectors  in 
space . 

There  are  currently  two  primary  methods  for  cooling 
satellite  systems  in  space:  (1)  evaporation  of  expendable 
cryogens  (liquid  helium,  nitrogen,  ammonia-methane,  and 
solid  hydrogen,  etc.),  and  (2)  closed  cycle  refrigerators 
(Stirling  cycle,  Vuilleumier  cycle,  etc.),  which  involve 
large    reciprocating    piston   masses   operating    at    low 
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frequency.  [Walker  (1983)  and  Smith,  et  al.  (1984)]  The 
disadvantages  of  these  current  cooling  methods  are  their 
short  lifetimes  (expendable  cryogens)  and  high  vibration 
levels  and  low  reliability  (closed  cycle  refrigerators).  The 
advantages  of  the  thermoacoustic  refrigerator  over  these  and 
other  systems  lies  in  its  simplicity  and  reliability.  The 
STAR  has  no  sliding  seals,  efficient  performance,  low  vibra- 
tion levels,  and  (hopefully)  a  longer  life  span. 

The  thermoacoustic  refrigerator  needs  to  be  surroun- 
ded by  a  vacuum  and  insulation  material  to  minimize  heat 
loss  due  to  thermal  conduction  and  radiation.  The  vacuum  of 
space  will  provide  the  best  insulation  to  determine  the 
refrigerator's  ultimate  efficiency.  The  absence  of  gravity 
removes  the  possibility  of  thermal  convection  as  an  addi- 
tional nuisance  heat  transport  mechanism. 
2.   Get  Away  Special  (GAS)  Program 

NASA's  GAS  program  allows  small,  self-contained 
payloads  to  fly  on  the  Space  Shuttle  in  GAS  canisters  for 
relatively  low  cost  ($10,000)  (Get  Away  Special  Team,  1984). 
The  GAS  can  is  five  cubic  feet  in  volume  and  can  house  a 
payload  of  up  to  200  pounds.  Each  payload  must  contain  its 
own  electrical  power,  control,  data  acquisition  and  storage 
facilities,  etc.  The  Shuttle's  astronauts  will  turn  a  switch 
on  or  off  at  designated  times  during  the  flight,  but  are 
otherwise  not  involved  with  the  experiment. 
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3.   Shared  Subsystems  from  NASA  G-313 

Another  group  of  NPS  students  and  faculty  (Boyd, 
et  al.,  1987)  have  taken  advantage  of  the  GAS  program  to 
measure  the  resonant  acoustic  modes  of  the  shuttle  payload 
bay  and  the  ambient  acoustic  environment  produced  as  a 
result  of  main  engine  and  booster  operation  during  launch. 
This  experiment  is  titled  "The  Space  Shuttle  Cargo  Bay 
Vibroacoustics  Experiment"  and  is  designated  by  NASA  as 
payload  G-313.  Several  subsystems  that  were  developed  for 
NASA  G-313  will  be  used  by  NASA  G-337  (STAR).  A  schematic  of 
the  STAR  in  its  GAS  can  is  shown  in  Figure  1-4.  One  of  the 
systems  borrowed  from  NASA  G-313  is  the  electronics  system 
used  to  run  the  experiment  and  record  the  data.  The  recorder 
system  consists  of  the  INTEL  model  BPK  5V75  magnetic  bubble 
memory  module  and  an  NSC  800  microprocessor-based  control- 
ler. Two  other  NPS  students,  LT  Charles  B.  Cameron,  USN  and 
CPT  Ronald  Byrnes,  USA,  will  be  designing  the  analog  elec- 
tronics and  software  to  integrate  these  systems  into  the 
STAR  experiment  for  their  master's  thesis.  A  block  diagram 
of  the  electronics  is  shown  in  Figure  1-5. 

Another  borrowed  system  is  the  power  supply,  which 
consists  of  Gates  brand  lead  acid  battery  cells  {five 
ampere-hour,  two  volts  each).  These  gelled  electrolyte  bat- 
teries are  ideal  for  the  STAR  due  to  their  high  power  den- 
sity, low  cost,  and  the  absence  of  outgassing  during  dis- 
charge  cycles.   NASA  G-313   used  a  one  layer  battery  of   68 
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cells  providing  680  watt-hours  of  energy  and  weighing  about 
80  pounds  (including  the  cells'  support  structure).  We  will 
be  using  two  battery  layers  with  as  many  cells  as  we  can,  up 
to  136.  This  will  give  us  a  minimum  of  680  watt-hours  and  a 
maximum  of  1,360  watt-hours  of  available  electrical  energy, 
and  a  battery  weight  between  80  and  160  pounds.  The  number 
of  battery  cells  we  can  use  will  depend  on  the  total  weight 
of  the  driver-resonator  assembly  and  its  auxiliary  equipment 
(vacuum  can,  gas  reservoir,  etc.)  and  electronics. 

The  use  of  the  GAS  can  imposes  certain  restrictions. 
Since  we  are  using  batteries  to  supply  the  power,  the  refri- 
gerator has  to  be  energy  efficient.  Also,  the  GAS  can  setup 
requires  the  STAR  to  be  compact  and  lightweight.  Figure  1-6 
shows  a  photograph  of  Hofler's  prototype  refrigerator.  This 
setup  is  approximately  six  feet  high.  In  comparison,  the 
maximum  payload  height  for  the  GAS  can  is  28.25  inches,  or 
less  than  2.5  feet.  These  considerations  played  a  major  role 
in  the  choice  of  equipment  for  and  the  design  of  the  STAR. 

C .   SCOPE 

Chapter  II  discusses  basic  acoustical  and  thermodynamic 
theory  as  applied  to  STAR.  The  results  of  this  theory  deter- 
mine the  driver  requirements.  The  theory  of  the  selected 
electrodynamic  driver  is  given  next.  Following  the  driver 
theory  is  a  presentation  and  discussion  of  a  computer 
program   that  models  the  coupled  system,   which  consists   of 
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Figure  1-6.   Photograph  of  Hofler's  (1986) 
Prototype  Refrigerator. 
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the   electrodynamic  driver  and  the  resonator  containing   the 
gas  mixture  and  stack . 

Chapter  III  describes  the  technique  for  measuring  the 
electrodynamic  driver  parameters.  These  parameters  are  used 
as  input  for  the  driver  design,  the  computer  model,  and  the 
efficiency  calculations.  This  measurement  technique  was 
developed  using  an  Altec  290-16K  electrodynamic  driver  and 
was  applied  to  the  final  driver. 

Chapter  IV  describes  the  measurement  of  the  STAR  driver 
parameters.  The  driver  chosen  for  STAR  is  a  custom  unit, 
built  by  Harmon-JBL!  <"^»^>  ,  which  uses  neodimium-iron-boron 
(NIB)  magnets.  The  STAR  driver  is  a  minor  modification  of  a 
new  line  of  compression  drivers  using  NIB  magnets  and  has 
characteristics  similar  to  the  JBL  2450 J . 

Chapter  V  discusses  the  construction  of  the  STAR  driver, 
including  the  housing  and  accessories  (reducer,  stiffener, 
bellows,  accelerometer ,  microphone,  electrical  feed- 
throughs ,  etc.)  necessary  to  integrate  it  with  the  resonator 
and  the  GAS  cannister  system.  Dimensional  shop  drawings  for 
the  parts  are  included  in  Appendix  C. 

Chapter  VI  gives  conclusions  and  recommendations  for 
further  development. 
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II.   THEORY 

A.   THERMOACOUSTIC  THEORY 

Thermoacoustic  theory  has  been  developed  in  detail  by 
Rott  (1969,  1974,  1975  and  1980)  and  adapted  to  the  thermo- 
acoustic refrigerator  by  Wheatley  et  al .  (1982,  1983a, 
1983b,  1935  and  1986),  Wheatley  and  Cox  (1985),  Hofler 
(1986)  and  Swift  (1989). 

The  space  thermoacoustic  refrigerator  (STAR)  basically 
consists  of  an  acoustic  driver  producing  sound  waves  in  a 
resonant  tube  (see  Figure  1-2).  This  resonator  is  filled 
with  a  mixture  of  helium-xenon  gas  (12.5%  xenon)  pressurized 
to  ten  atmospheres ,  and  contains  a  stack  of  plastic  plates . 
This  Chapter  will  present  a  qualitative  model  for  the  ther- 
moacoustic heat  pumping  process  followed  by  a  quantitative 
development  for  heat  and  work  flow  at  a  plate.  The  Chapter 
concludes  with  the  efficiency  of  a  stack  of  plates  that  are 
very  much  shorter  than  one  quarter  of  a  wavelength. 

1 .   Qualitative  Picture 

Consider  a  parcel  of  gas  that  moves  back  and  forth 
along  one  of  the  plates  at  the  acoustic  frequency  (see 
Figure  II-l).  As  it  moves,  the  gas  parcel  will  experience 
changes  in  temperature  and  volume.  Part  of  the  temperature 
changes  come  from  the  adiabatic  compression  and  expansion  of 
the  gas  by  the  sound  pressure,   and  part  as  a  consequence  of 
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the  local  temperature  of  the  plate  itself.  A  temperature 
gradient  may  develop  along  the  plate  as  a  result  of  the 
operation  of  the  refrigerator.  The  temperature  and  volume 
changes  can  be  described  by  six  separate  steps  (the  changes 
in  volume  are  indicated  by  the  change  in  size  of  the  square 
parcel  in  the  figure) . 

Assume  the  plate  is  at  a  uniform  temperature,  T.  The 
parcel  starts  at  position  1  (X  =  0)  with  a  temperature  of  T. 
The  sound  wave  moves  the  parcel  to  the  right  to  position  2. 
The  parcel  has  undergone  adiabatic  compression  and  its  tem- 
perature is  now  T"*"*"  .  Since  the  temperature  of  the  gas  is 
higher  than  that  of  the  plate,  heat  (O;  will  flow  from  the 
gas  parcel  to  the  plate.  The  parcel's  volume  is  decreased 
and  its  temperature  lowered  to  T"*" .  This  heat  flow  also 
causes  the  plate's  temperature  to  increase  at  the  position 
X"^ .  Parcel  position  3  is  actually  in  the  same  location  as 
position  2,  but  it  is  displaced  vertically  in  the  figure  for 
clarity.  The  parcel  now  moves  to  position  4  and  expands 
adiabatically  to  a  new  temperature  of  T-  ^  repeat  of  this 
process  puts  the  parcel  in  position  5  with  a  temperature  of 
T  .  Now  the  temperature  of  the  gas  is  lower  than  that  of 
the  plate  and  heat  (Q')  flows  from  the  plate  to  the  gas 
parcel,  expanding  the  parcel  and  raising  its  temperature  to 
T~~.  This  flow  causes  the  plate's  temperature  to  decrease  at 
position  X- .  The  parcel  moves  to  the  right  again  under 
adiabatic   compression   and  we  are  back  where  we  started   in 


position  1  with  temperature  T  (adapted  from  Wheatley,  et 
al. ,  1985) . 

This  gas  parcel  cycle  is  repeated  at  the  resonator's 
operating  frequency.  For  the  STAR  this  frequency  is  between 
approximately  250  and  600  cycles  per  second,  depending  on 
the  gas  mixture  (specifically  the  speed  of  sound  in  the  gas 
mixture)  and  the  length  and  shape  of  the  resonator. 

Notice  that  the  temperature  of  the  gas  parcel  at 
position  X  =  0  is  different  depending  on  which  direction  the 
gas  parcel  is  moving.  It  ±s  this  phase  shift  in  temperature 
relative  to  motion  that  produces  the  thermoacoustic  effect, 
as  we  will  show  in  the  calculations  that  follow. 

It  is  also  important  to  note  that  in  order  for  heat 
to  flow  between  the  gas  parcel  and  the  plate,  the  parcel 
must  be  vertically  located  within  about  a  thermal  penetra- 
tion depth  of  the  plate.  The  thermal  penetration  depth  (Sk) 
is  the  distance  that  heat  can  diffuse  through  the  fluid 
during  a  time  1/a^,  where  <u  is  the  acoustic  angular 
frequency . 

Wheatley  et  al .  (1986)  describe  thermoacoustic  en- 
gines as  consisting  of  long  trains  of  these  gas  parcels,  all 
about  a  thermal  penetration  depth  from  the  plate.  The  par- 
cels draw  heat  from  the  plate  at  one  extreme  of  their  oscil- 
latory motion  and  deposit  heat  at  the  other  extreme.  Adja- 
cent heat  flows  cancel  except  at  the  ends.  The  net  result  is 
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that  an  amount  of  heat  Q  is  passed  from  one  end  of  the  plate 
to  the  other. 

2.   Single  Plate 

Swift  (1989)  illustrated  the  basic  principles  of  the 
thermoacoustic  engine  by  examining  a  simple  example  of  a 
single  plate  in  a  gas-filled  tube  such  as  we  described 
above.  In  this  example  the  acoustic  and  thermodynamic  ef- 
fects are  nearly  distinct.  In  the  absence  of  the  plastic 
plate,  the  plane  standing  wave  sustained  in  the  tube  has 
well-known  acoustic  properties.  It's  pressure  and  x-direc- 
tion  velocity  (along  the  length  of  the  tube)  are  given  as 
(using  peak  values): 

Pi  =  p«  sin(x/A)  ( II-l ; 

and      Uj  =  i(p«/p„a)  cos(x/A)  (II-2) 

where  P^  is  the  pressure  amplitude  at  the  pressure  antinode, 
the  reduced  wavelength  A  =  ^/27T  =  a/a>,  ?\  is  the  wavelength, 
a  is  the  speed  of  sound  in  the  gas,  <m  is  the  angular  fre- 
quency, and  p„  is  the  mean  density  of  the  gas.  The  subscript 
1  indicates  the  first  order  time  oscillation,  and  the  sub- 
script m  indicates  the  mean  value.  The  total  pressure  can  be 
given  as 

P  =  Pm  +  pie^*^*-.  (II-3) 

A  similar  expression  represents  the  total  velocity.  Follow- 
ing Swift,  we  will  assume  that  the  first  order  in  the  acous- 
tic  amplitude   is  sufficient  for  all  acoustic   and   thermo- 
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dynamic  variables  (pressure,  velocity,  temperature,  density, 
and  entropy) . 

The   sound  wave  is  adiabatic  and  has  an   oscillatory 
temperature  given  by: 

Tx    =    (T„e/PmC,-)  px,  (11-4) 

where  3  =  -  (dp/dT)F=/Pm  is  the  isobaric  thermal  expansion 
coefficient  and  Cp.  is  the  isobaric  (constant  pressure)  heat 
capacity  per  unit  mass.  Note  that  T^  and  pi  are  in  phase. 
For  ideal  gases,  thermodynamics  tells  us  that 

Tme/PmC,=.  =  (^-l)Tm/;^p„,  (II-5) 

where  -i    =  the  ratio  of  isobaric  to  isochoric  specific  heats 

(5/3  for  monatomic  gases,  smaller,  but  greater  than  one,  for 

other   gases).   Combining  Equations  II-4  and  II-5  gives   us : 

Ti/T„  =  [{-6-l)/-6]    pi/p™.  (II-6) 

The   "'.ntroduction   of   the  plastic   plate   into   the 

standing  wave  modifies  the  original,  unperturbed  temperature 

oscillations.   This   modification   is  due  to  the   heat   flow 

between  the  gas  and  the  plate,   as  described  previously.  The 

temperature   is  modified  in  both  magnitude  and  phase  for  gas 

about   a   thermal   penetration  depth  away   from   the   plate. 

According  to  Swift,   this  results  in  two  important   effects: 

1)   a  time-aveiage  heat  flow  near  the  surface  of  the   plate, 

along   the   direction   of  acoustic   vibration,   and   2)   the 

generation   or   absorption  of  real  acoustic  power   near   the 

surface   of   the  plate.   In  refrigerators  the  effect   is   an 

absorption   of  acoustic  power.   In  prime  movers  the  acoustic 
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power  is  generated  by  the  temperature  gradient  present  in 
the  plate. 

There  are  several  assumptions  we  will  make  to 
simplify  calculations,  again  following  Swift.  We  assume  that 
the  plate  is  short  enough  compared  to  a  reduced  wavelength 
(Ax  «  A)  and  far  enough  from  both  velocity  and  pressure 
nodes  that  p,^  and  Uj.  can  be  considered  uniform  over  the 
entire  plate.  We  assume  the  gas  has  zero  viscosity,  so  that 
Ui  does  not  depend  on  y  (we  had  already  made  this  assumption 
by  looking  at  Uj.  in  the  x-direction  only)  .  We  assume  that 
the  plate  has  a  large  enough  heat  capacity  per  unit  area 
that  its  temperature  does  not  change  appreciably  at  the 
acoustic  frequency.  We  assume  that  the  plate  has  a  mean 
temperature  gradient  in  the  x-direction  vTm.  Finally,  we 
neglect  the  plate's  and  gas'  thermal  conductivity  in  the  x- 
direction . 

Applying   these   assumptions ,   we  see  that  the   mean 

fluid  temperature  (T,„(x))  is  the  same  as  that  of  the  plate. 
Swift  calculated  the  oscillating  fluid  temperature  using 
the  general  equation  of  heat  transfer.  He  kept  only  first- 
order  terms,  neglected  thermal  conduction  along  x,  and 
applied  the  boundary  condition  Ti(0)=0  imposed  by  the  plate. 
The  resulting  equation  is: 

Ti  (y)  =  [  (T,„a/p„c^)Pi-OT„/.:^)u^]  [i_exp[-{  l  +  i)y/S]  ]  .  (II-7) 
To  interpret  this  equation,   we  look  at  it  in  the  limit  that 
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the  gas  is  far  enough  from  the  plate  (y»SK)  to  make  negli- 
gible thermal  contact  with  the  plate.  This  gives: 

Ti(y)-i^[  (T„e/p„c^)pa-(vT™/ai)Ui]  (II-3) 

The    first   term   in   Equation   II-8   is   the    temperature 

oscillation  due  to  the  adiabatic  compressions  and  expansions 

of  the  gas  (see  Equation  II-4) .   The  second  term  comes   from 

the   mean   temperature   gradient   in  the   gas .   As   the   gas 

oscillates  along  x  with  displacement  amplitude  ui/'jj,  the 
temperature   at   a   given  point  in  space   oscillates   by   an 

amount  sJ-Tm  Ui/-x-  even  if  the  temperature  oi  a  given  piece  of 
fluid  remains  constant.  The  actual  temperature  oscillations 
are  a  linear  superposition  of  these  two  effects. 

The  y  dependent  part  of  equation  II -7  is  complex.  It 
approaches  1  for  y  »  Sk,  as  given  previously.  It  approaches 
zero  for  y  «  Sk ,  where  the  plate  imposes  the  condition 
Ti  =  0.  Most  importantly,  for  y  5:  Sk  ,  its  magnitude  is  still 
of  the  order  of  1,  but  it  has  a  substantial  imaginary  part. 
This  phase  shift  in  the  oscillating  temperature  of  the 
standing  wave  at  y  a;  Sk  ,  due  to  the  thermal  presence  of  the 
plate,  is  an  important  result  because  it  leads  directly  to 
the  time-averaged  heat  flow  in  the  x-direction.  It  is  the 
same  phase  shift  we  emphasized  in  the  description  of  the  gas 
parcel  motion. 

a.   Heat  Flow 

Swift  argues  that  since  we  are  neglecting   ordi- 
nary  thermal  conductivity  in  the  x  direction,   the  only  way 
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heat  can  be  transported  along  x  is  by  the  hydrodynamic 
transport   of  entropy,   carried  by  the  oscillatory   velocity 

q=   =  Tm  Pm  STTTT  =  1/2  PmC.»  Im  [Tx]  Ux.       (II-9) 
The   line  above  the  quantity  SxUx  indicates  the  time-average 

of  the  first-order  entropy  and  velocity  product.  The  heat 
flow  is  a  second  order  quantity,  signified  by  the  subscript 
2. 

The  total  heat  flow  C>~,   along  the  plate,   in  the  x- 
direction,  is  found  by  integrating  q^  over  the  y-z  plane: 

®=  =  TT  oX*^q=  dy  (11-10) 

where   IT   is  the  perimeter  of  the  plate  in   the   y-z   plane. 

Substituting  for  Tj_  and  performing  the  integration  gives: 

Qz  =  -1/4  rT3K{T„,e)PxUx(r-l)  (11-11) 

where  ITSk  is  the  thermodynamically  active  area  in  a  plane 
perpendicular  to  the  longitudinal  acoustic  motion,  Tm0  is 
the  heat  parameter  of  the  gas,  and  T  is  defined  as  the  ratio 
of   actual  temperature  gradient  to  the   critial   temperature 

gradient  (^T/vT<=r-ifc)  .  The  critical  temperature  gradient 
occurs  when  the  temperature  change  along  the  plate  just 
matches  the  temperature  change  due  to  the  adiabatic 
compression  of  the  gas,  and  no  heat  flows  between  the  gas 
and  the  plate.  It  is  the  boundary  between  the  refrigerator 
and  prime  mover  functions  of  the  thermoacoustic  engines. 
Equation  11-11  shows  that  when  f  <  i,  heat  flows 
up   the   temperature   gradient  from  cold   to   hot   and   work 
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(acoustic  power)  is  absorbed,   as  for  a   refrigerator.   When 
r  =  1  there  is  no  heat  flow.  When  T  >  1,  heat  flows  down  the 

temperature   gradient  from  hot  to  cold  and  acoustic  power  is 
produced,  as  for  a  prime  mover. 

Note  that  the  total  heat  flow  is  proportional  to 
the  area  TTSk  ,  and  to  Tm3  (=1  for  ideal  gases)  .  It  is  also 
proportional  to  the  product  piUi,  and  so  equals  zero  if  the 
plate  is  at  either  a  pressure  node  or  a  velocity  node  of  the 
standing  wave.  The  maximum  value  of  the  product  occurs 
hal"^ay  between  the  nodes.  Finally,  the  heat  flow  is 
proportional  to  the  temperature  gradient  factor  T-i .  For 
"^T^  ^T^i-i*:. ,  F-l  0  and  the  heat  flow  is  toward  the  pres- 
sure node,  while  for  ^T„  ■■  ^T^r-±^-,  F-l  ■■:  0  and  the  heat  flow 
is  away  from  the  pressure  node.  If  suitable  heat  exchangers 
at  temperatures  Tm  and  Tc  are  installed  at  the  ends  of  the 
plate  (with  T„  -  t,^  =  vT„,  Ax),  this  heat  flow  carries  heat 
from  one  exchanger  to  the  other. 

The  heat  flow  is  small  under  ordinary 
circumstances.  However,  in  closed  resonators  it  is  possible 
to  achieve  sound  amplitudes  many  orders  of  magnitude  higher 
than    those    of   ordinary   conversation.    Since    O^    is 

proportional  to  p«^ ,  and  since  iii  practical  acoustic  engines 
the  entire  cross  section  of  the  standing  wave  is  filled  with 
plates  (spaced  roughly  4Sk  apart),  very  high  heat  flows 
and/or  a  large  T^^  -  t^  may  be  achieved. 


b.   Work  Flow 

The  work  flow  {i.e.   acoustic  power)  is  given  by 

the  work  per  cycle  times  the  rate  at  which  that  work   occurs 

(the  acoustic  frequency  f).  From  thermodynamics,  the  average 

acoustic  power  produced  per  unit  volume  is: 

^2  =  -im/p^)    ipiPi  =  -  1/2  u>0    pi  Im[Tx].    (11-12) 
The   gas  about  a  thermal  penetration  depth   away 

from  the  plate  "breathes,"  because  of  thermal  expansion  and 

contraction,   with   the   right   time  phase  with   respect  to 

oscillating   pressure  to  do  (or  absorb)  net   work.   This  is 

exactly   the   same  gas  that  we  have  seen  is  responsible  for 

the   heat  flow.   Gas  elsewhere  is  ineffective  in   doing  (or 

absorbing)  work.  The  density  oscillations  for  y  .^    ^^..    ^^^  ^-^ 

y   »   Sk  are  in  phase  with  the   pressure   oscillations,   and 

hence  do  (or  absorb)  no  net  work. 

The  total  acoustic  power  1^2  produced  is  found  by 
integrating  i^z  over  all  space,  as  with  the  heat  flow: 

Wz  =  1/4  TT-Sk  (-^-1)  (px  =  /?,„a)  (T-l)  Ax/A,    (11-13) 
where  all  the  terms  have  been  defined  previously. 

The  acoustic  power  is  proportional  to  the  volume 
ITSkAx  of  fluid  that  is  about  a  thermal  penetration  depth 
from  the  plate.  It  is  proportional  to  Pi^,  and  so  is  quad- 
ratic in  the  acoustic  amplitude  (as  was  the  heat  flow)  and 
vanishes  at  pressure  nodes.   Finally,   W_.  j_^    proportional  to 

(f-i),   the   same  temperature  gradient  factor  as  appeared  in 
O2  (Equation  11-11).  When  vT™  =  vT^^i,^,  (f-l)  =  0,  and  there 
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are  no  temperature  oscillations  in  the  fluid  other  than 
those  due  to  adiabatic  processes,   and  no  acoustic  power  is 

absorbed  or  generated.  For  vTm  >  vT^r-it^,  r  -i  >  0  and  acous- 
tic power  is  produced  near  the  plate.  Whether  this  power 
increases  the  amplitude  of  the  standing  wave,  is  radiated 
away  to  infinity,  is  simply  absorbed,  or  flows  through  an 
acoustic-to-electric  transducer  to  generate  electric  power 
depends  on  details  of  the  resonator,  not  on  the  plate  itself 

or  on  the  standing  wave  near  the  plate.  For  vTm  <  vT=r-i«-, ,  T- 
1    <  0  and  acoustic  power  is  absorbed  near  the  plate.   For   a 

tube   without   plates   which  has  a  diameter   less   than   the 

wavelength,  at  constant  temperature  ( vT  =  0),  this  work  flow 

is   responsible   for  the  ordinary  thermal  attenuation   of   a 

sound  wave  (Kinsler,  et  al . ,  1932). 

c.   Efficiency 

We  combine  Equations  11-11  and  11-13  to  get   the 

efficiency   of   the   plate  with  no  viscous   or   longitudinal 

conduction  losses: 

^  =  W2/O2  =  (v-i)/(T„3)  (<D  Ax  pi)/(Fma==Ui)  .  (11-14) 
Since  u^  =  Uo  sin  x/A,  and  pi  =  p„a  Uo  cos  x/X, 

n  =  {;i-l)/(T„e)  (Ax/A  tan  x/a)  .  (11-15) 

For  X  «  ?v,   n  =  (5-l)/{T„a)  (Ax/x),  (11-16) 

and      n  =  ncr-^^t/r  (11-17) 

This  efficiency  approaches  the  Carnot  efficiency 

as  the  power  output  and  heat  transfer  rates  approach  zero. 
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We  can  make  a  similar  calculation  for  the 
refrigerator  mode  of  operation,  where  the  relevant 
efficiency  is  the  coefficient  of  performance,  COP  =  O^/l^^. 
We  find  that 

COP  =  r  COP<=.^„«t.  (11-13) 

3.   Short  Stack 

After  developing  the  basic  principles  of  the  thermo- 
acoustic  engine  using  a  simplified  example,  Swift  developed 
a  more  realistic  model  by  including  viscosity,  longitudinal 
thermal  conductivity,  finite  (instead  of  infinite)  plate 
heat  capacity,  and  many  plates.  He  still  made  the  assump- 
tions that  the  temperature  spanned  is  much  less  than  the 
absolute  temperature.  that  the  length  of  the  plates  is  much 
shorter  than  the  reduced  wavelength,  and  that  the  spacing 
between  the  plates  is  greater  than  a  few  penetration  depths. 

Using  the  equation  of  motion  for  the  gas,  the  boun- 
dary condition  that  the  velocity  at  the  gas-plate  interface 
is  zero,  the  continuity  equation  of  the  gas,  the  heat  flows 
in  the  gas  and  plates,  and  basic  thermodynamic  properties, 
Swift  derived  a  wave  equation  for  pressure  Pi(x)  in  terms  of 
dTm/dx ,  material  properties  and  geometry.  Swift  also  derived 
an  equation  for  enthalpy  flow  along  x  in  terms  of  Pi(x), 
T„(x),  geometry,  and  material  properties  (including  Prandtl 
number  of  the  gas  cr  =  Cf.u/k,  where  u  is  the  viscosity,  and  k 
is  the  thermal  conductivity  of  the  gas). 
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In  the  boundary  layer  approximation,  for  a  short 
(Ax  «  ^)  stack  of  plates,  and  neglecting  viscosity,  the  heat 
flux  and  acoustic  power  are  given  by  expressions  very  simi- 
lar to  the  simple  results  of  the  single  plate  example. 
Practical  engines  can  be  expected  to  have  all  available 
cross-sectional  area  filled  with  plates  spaced  about  4.3k 
apart.  The  only  effect  of  longitudinal  (i.e.,  along  x) 
thermal  conductivity  is  to  add  to  the  heat  flow  in  the  x 
direction  by  simple  conduction.  For  gases  this  effect  is 
normally  negligible. 

The  inclusion  of  viscosity  adds  considerable 
algebraic  complication,  and  a  little  conceptual  difficulty 
as  novj  ui  is  a  function  of  y. 

Without  viscosity,  vTcr-n=.  is  the  temperature 
gradient  for  which  the  temperature  oscillation  at  a  given  x- 
location  is  zero.  But  in  the  presence  of  viscosity,  T^ 
depends  on  x  as  well  as  y,  so  no  well-defined  vTc-it^  exists. 
To  the  lowest  order  in  -fU  (where  a  is  the  fluid's 
Prandtl  number),   the  hydrodynamic  heat  flow  is  increased  by 

l/ll-Vo").  This  factor  arises  because  the  mean  velocity  <Ui> 
underestimates  the  velocity  with  which  entropy  is  convected 
along  the  stack.  As  we  saw  in  the  single  plate  calculations, 
the  convective  entropy  transport  occurs  mostly  at  a  distance 
Sk  from  the  plate.  The  velocity  there  is  higher  than  the 
mean  velocity  by  1/(1-t/o'). 
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The  lowest  order  effect  of  viscosity  on  the  acoustic 
power  is  the  power  dissipated  by  viscous  shear  in  the  fluid 
within  the  boundary  layer,  a  well-known  fluid  mechanical 
result . 

For  accurate  calculations,  the  full  general  equation 
must  be  used;  but  for  rough  estimates,  the  results  presented 
here  are  good  approximations,  and  are  much  simpler  to 
compute  with.  In  addition,  the  results  presented  here  are 
fairly  easy  to  interpret  physically.  The  expressions  for 
arbitrary  viscosity  are  complicated,  but  those  to  lowest 
order  in  viscosity  are  simple.  It  is  easy  to  see  the  single 
plate  expressions  for  heat  flow  and  acoustic  power  about  a 
thermal  penetration  depth  from  the  plates,  modified  as 
longitudinal  thermal  conductance  adds  to  the  heat  flow,  as 
viscous  shear  dissipates  acoustic  power,  and  as  non- ideal 
plate  properties  modify  the  thermal  boundary  condition  at 
the  fluid-plate  interface. 

4 .   Design  Considerations 

One  of  the  important  considerations  in  the  design  of 
a  thermoacoustic  refrigerator  is  the  location  of  the  stack 
of  plates  with  respect  to  the  pressure  and  velocity  nodes 
and  antinodes  of  the  acoustic  standing  wave  in  the  resonant 
tube.  Experiments  performed  by  Wheatley,  et  al.  {1986) 
showed  that  the  heat  flow  as  a  function  of  position  fit  a 
simple  sine  curve  whose  spacial  period  is  half  the  wave- 
length of  the  acoustic  standing  wave.  By  noting  how  the  sign 


33 


of  the  temperature  difference  varied  with  respect  to  the 
plate's  location  in  the  sound  wave.  they  saw  that  heat 
always  flows  in  the  direction  of  the  closest  pressure  anti- 
node.  This  effect  is  expected  because  a  parcel  of  gas  moving 
in  the  direction  of  a  pressure  antinode  is  compressionally 
warmed  and  will  transfer  heat  to  the  plate.  A  parcel  moving 
toward  a  pressure  node  is  cooled  by  expansion  and  will  draw 
heat  from  the  plate.  At  both  the  pressure  nodes  (velocity 
antinodes)  and  antinodes  (velocity  nodes)  heat  flow  drops  to 
zero.  Thus  the  acoustic  heat  flow  depends  on  both  the  acous- 
tic pressure  and  the  fluid  velocity. 

As  a  plate  or  stack  of  plates  is  moved  away  from  the 
pressure  antinode,  the  temperature  gradient  developed 
becomes  smaller.  At  a  quarter  of  a  wavelength,  no  gradient 
forms  (the  critical  temperature  gradient  equals  zero).  This 
positioning  effect  is  important  in  the  design  of  a  refriger- 
ator, because,  together  with  the  length  of  the  plates.  it 
places  an  upper  limit  on  the  maximum  temperature  drop 
possible  across  the  stack. 

B.   DRIVER  REQUIREMENTS 

In  order  to  drive  the  thermoacoustic  refrigerator  we 
needed  a  device.  a  transducer,  that  would  convert  the  elec- 
trical energy  of  a  battery  pack  to  acoustic  energy,  as 
discussed  in  Subsection  I. A. 3.  We  wanted  as  high  a  pressure 
(or   force)  and  volume  velocity  (or  motion)  as   possible   to 
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get  the  most  acoustic  power  provided  to  the  system.  There 
were  several  types  of  electroacoustic  transducers  to  choose 
from:  electro-static,  electret,  electrodynamic ,  magneto- 
strictive,  electro-strictive ,  and  piezoelectric.  Only  one  of 
these,  the  electrodynamic  transducer,  can  provide  our  system 
with  both  the  large  forces  and  velocities  necessary.  This 
transducer  will  be  described  in  the  next  section.  The 
remainder  of  this  section  describes  the  rest  of  the  trans- 
ducers and  explains  why  each  was  not  appropriate  for  use  in 
the  thermoacoustic  refrigerator. 

The  electrostatic  transducer  consists  of  a  pair  of 
charged  electrodes,  or  capacitor  plates,  one  of  which  is 
held  stationary  while  the  other,  the  diaphragm,  moves  in 
response  to  electrical  excitation.  The  electret  is  a  kind  of 
electrostatic  transducer  that  has  a  diaphragm  of  polarized 
plastic.  and  therefore  does  not  need  an  external  voltage 
supply  to  the  diaphragm  (Kinsler,  et  al,  1932,  p.  350).  If 
an  alternating  voltage  is  applied  across  the  plates,  the 
diaphragm  moves,  thereby  radiating  an  acoustic  wave. 
Unfortunately  there  is  very  little  force  (or  pressure)  and 
displacement  associated  with  this  motion.  This  type  of 
transducer  may  be  useful  as  a  microphone  but  is  not  an 
appropriate  choice  to  drive  the  thermoacoustic  refrigerator. 

In  a  magnetostrictive  transducer  a  change  in  the 
magnetic  polarization  in  the  material  causes  an  elastic 
strain  (Wilson,   1985,   p.  3.).  Although  the  forces  are  more 
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than  adequate  the  strain  produces  too  little  motion  to  be 
useful  to  us. 

In  electrostrictive  and  piezoelectric  transducers  an 
externally  applied  electric  field  causes  a  change  in  the 
dielectric  polarization  of  the  material,  which  in  turn 
causes  an  elastic  strain  (Wilson,  1985,  p.  3).  Once  again 
there  is  sufficient  force  but  the  displacement  due  to  the 
strain  is  too  small  to  be  useful . 

Because  these  techniques  have  sufficient  force,  it  would 
be  possible  to  overcome  the  displacement  limitation  using 
mechanical  leverage  or  other  forms  of  mechanical  impedance 
transformation,  but  this  was  not  considered  due  to  the 
associated  increase  in  device  complexity. 

C.   ELECTRODYNAMIC  THEORY 
1 .   Physics 

The  electrodynamic  driver  is  a  moving  coil 
transducer  which  converts  electrical  energy  to  kinetic  ener- 
gy (motion) .  It  consists  of  a  diaphragm  attached  to  a  cylin- 
drical coil  of  wire  (voice  coil)  that  is  suspended  in  a 
fixed  magnetic  field  (see  Figure  II-2).  If  an  alternating 
current  is  supplied  to  the  coil,  the  interaction  of  the 
current  and  the  magnetic  field  will  induce  a  force  on  the 
coil  so  that  the  diaphragm  moves.  The  magnitude  of  the  force 
is  F  =  BlI,  where  B  is  the  magnitude  of  the  magnetic  induc- 
tion,  1  is  the  length  of  the  coil  and  I  is  the  current.  The 
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Figure  II-2. 


Schematic  of  an  Electrodynamic  Transducer 

Showing  (A)  Magnet,  (B)  Back  Plate, 

(C)  Pole  Piece,  (D)  Acoustic  Seal, 

(E)  Diaphragm,  (F)  Voice  Coil, 

(G)  Surround  (Suspension) ,  and 

(H)  Electrical  Lead  Connection. 
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direction  of  the  force.  given  by  the  cross  product  of  the 
current  vector  and  the  magnetic  field  (from  Ampere's  law  F  = 

II  X  B) ,  is  everywhere  orthogonal  to  the  current  and  the 
radial  magnetic  field.  The  motion  of  the  diaphragm  produces 
sound  waves  with  a  frequency  equal  to  that  of  the  alterna- 
ting current.  Another  important  property  of  the  electrodyna- 
mic  transducer  is  that  the  motion  of  the  voice  coil  induces 
an  electromotive  force  (emf)  which  equals  Bl  times  the 
velocity  of  the  voice  coil  (emf  =  Blv).  This  "back  emf" 
increases  the  electrical  impedance  of  the  driver  and  can  be 
used   to  monitor  the  driver  motion  as  discussed   in   Section 

I I I  .  B  . 

The  driver  is  attached  to  the  resonator  via  a 
reducer  and  bellows.  The  reducer  is  a  cone-shaped  piece  of 
aluminum  mechanically  connecting  the  driver  voice  coil  which 
produces  the  force  to  the  bellows .  which  in  turn  is  in 
contact  with  the  gas  in  the  resonator.  The  reducer  is 
designed  to  be  lightweight,  rigid,  and  strong  for  a  direct 
transfer  of  motion  and  energy  from  the  moving  coil  to  the 
bellows.  The  bellows  is  a  lightweight  and  flexible  gastight 
seal  between  the  driver  and  the  tube.  An  aluminum  plate  is 
used  as  an  interface  between  the  reducer  and  the  bellows  to 
keep  the  bellows  face  rigid.  For  the  STAR  driver  the  reducer 
and  plate  are  incorporated  into  a  single  unit  machined  from 
aluminum  bar  stock. 
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Figure  II-3  shows  a  mechanical  model  of  the 
electrodynamic  transducer  as  a  forced  simple  harmonic 
oscillator  (SHO) ,  which  consists  of  a  moving  mass  (m^ ) ^ 
stiffness  (k),  mechanical  damping  or  resistance  (Rm),  and 
driving  force  (F).  Included  in  the  moving  mass  are  the 
masses  of  the  voice  coil,  reducer,  stiffener,  and  bellows. 
The  stiffness  includes  the  suspension  stiffness  of  the  voice 
coil  surround,  the  stiffness  of  the  gas  volume  trapped 
behind  the  diaphragm,  and  the  stiffness  of  the  bellows.  The 
equations  governing  the  SHO  can  thus  be  used  to  determine 
the  moving  mass  and  stiffness  of  the  driver.  The  primary 
relationship  used  is: 

"^-^^   =  k/m^,  (11-19) 

or,   the   resonant   angular   frequency   squared   equals   the 

stiffness   divided  by  the  moving  mass.   The  use  of  this   and 
other   relationships   will  be  discussed  further   in   Chapter 

III,  Driver  Parameter  Measurements. 
2 .   Equivalent  Circuit 

An  electrical  circuit  can  be  drawn  that  is 
equivalent  to  an  electromechanical  system  if  the  values 
given  for  Rm,  Lm,  and  Cm  in  the  circuit  produce  the  same 
electrical  effects  as  the  electromechanical  system  itself. 
In  the  case  of  the  electrodynamic  driver,  the  equivalent 
circuit  is  shown  in  Figure  II-4  with  the  values  given  by 
Kinsler,  et  al .  (1982)  as: 
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Electrodynamic  Transducer--Forced  Simple 
Harmonic  Oscillator  Model. 
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Figure    II-4. 


Equivalent  Electrical  Circuit  for  the 
Electrodynamic  Transducer. 
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Rm  =  (Bl)=/R„,  (11-20) 

L„  =  (Bl)=/k,  (11-21) 

C„  =  in«/(Bl)  =  ,  (11-22) 

where   81   is   the   transduction   coefficient,    R„   is   the 

mechanical   resistance,   and  k  and  m^  are  the  stiffness   and 

moving  mass  as  described  previously. 

The  properties  of  the  equivalent  circuit  allow  us  to 
evaluate  the  performance  of  the  electrodynamic  driver  as  if 
it  were  an  entirely  electrical  system.  The  amount  of  power 
the  driver  generates  vm_11  depend  on  the  load  it  sees.  Since 
one  of  the  primary  purposes  of  this  thesis  is  to  maximize 
the  ratio  of  the  acoutic  power  delivered  to  the  electrical 
power  provided  (the  electroacoustic  efficiency),  it  is 
necessary  to  provide  a  load  to  the  system  in  order  to 
evaluate  driver  performance.  This  load  is  originally 
represented  as  a  complex  electrical  impedence  in  the 
equivalent  circuit.  The  actual  value  of  the  load  will  depend 
on  the  design  of  the  resonator  but  can  be  characterized  by 
an  acoustic  impedance  which  can  be  transformed  to  an 
equivalent  electrical  load  impedance. 
3 .   Losses 

In  our  search  for  the  highest  electroacoustic 
efficiency  it  is  useful  to  understand  where  possible  losses 
may  occur  and  do  what  can  be  done  to  minimize  them.  Losses 
that  are  an  intrinsic  part  of  the  electrical  driver  are  eddy 
current   losses   in  the  magnetic  structure  generated  by   the 


42 


current  in  the  voice  coil.  Joule  heating  losses  due  to  the 
drive  current  dissipated  by  the  voice  coil  electrical 
resistance,  and  mechanical  losses  due  to  the  internal  fric- 
tion of  the  surround  and  the  viscosity  of  the  gas  in  the 
"gap."  These  losses  need  to  be  taken  into  account.  Very 
little  can  be  done  about  the  eddy  current  losses  and  the 
viscosity  of  the  gas  which  is  necessary  to  improve  heat 
removal  from  the  voice  coil.  The  losses  due  to  surround 
(suspension)  mechanical  resistance  can  be  controlled  by 
choice  of  the  materials  and  removal  of  extraneous  suspension 
material  (see  Section  III.H). 

D.   COUPLED  SYSTEM 

1 .   Computer  Model 

Any  discussion  of  the  driver  section  of  the 
thermoacoustic  refrigerator  is  incomplete  without  some 
understanding  of  how  the  driver  parameters  are  coupled  to 
the  resonator  (acoustic  load).  The  equivalent  electro- 
mechanical circuit  for  the  coupled  system  is  the  same  as 
that  shown  in  Figure  II-4  with  the  acoustic  load  impedance: 

Zi  =  -Pccsj  [  (l  +  ja/k)/(l+(a/k)==)  ♦        (11-23) 
(  (cos(kL)sin(kL)  +  j  sinh  ( <xL)  cosh(«L)  )  ^ 

(sin=(kL)cosh==(<xL)  +  cos  =  (  kL)  sinh^(  <xL)  )  )  ]  , 
where  p^  is  the  density  of  the  gas,   c  is  the  speed  of  sound 

(labeled  "a"  in  previous  discussions),  s  is  the  cross- 
sectional   area   of   the  resonator,   k  is   the   wave   number 
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{=  tu/c)  ,  s  =  u}o/(2Qc)  ,  Q  is  the  quality  factor  of  the  reso- 
nator, and  L  is  the  effective  length  of  the  resonator 
(Kinsler,  et  al.,  1982).  This  circuit  combines  the  equiva- 
lent circuit  of  the  electrodynamic  transducer  with  the  com- 
plex impedance  load  representing  the  resonator.  One  way  to 
examine  this  coupling  is  through  computer  modeling. 

LT  Michael  P.   Susalla,   USN ,   developed  a   computer 
program   which  models  the  thermoacoustic  refrigerator  as  the 

equivalent  electrical  circuit.  The  program  models  each  of 
the  mechanical  and  acoustical  components  as  equivalent 
electrical  components,  as  discussed  above.  Using  the 
electrical  equation  from  the  combined  circuit  we  determine 
the  system  performance  as  different  parameters  are  varied, 
a.   Input 

There  are  two  sets  of  input  parameters  needed 
for  the  model:  (1)  driver  parameters  and  (2)  resonator 
parameters.  The  driver  electrical  parameters  include  the 
voice  coil  inductance  (LE)  and  resistance  (RE),  and  the 
current  (I)  supplied  to  the  voice  coil  (which  is  held  at  one 
ampere  for  all  test  cases).  The  driver  mechanical  parameters 
used  in  the  model  are  the  same  equivalent  resistance,  induc- 
tance, and  capacitance  as  those  developed  for  the  equivalent 
circuit.  Therefore  the  parameter  values  needed  for  input  are 
the  transduction  coefficient  ( BL) ,  the  mechanical  resistance 
(RM),   the  moving  mass  (M) ,  the  suspension  stiffness  (SSUP), 
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the   volume  of  the  gas  behind  the  diaphragm  (VOL) ,   and   the 
cross-sectional  area  of  this  volume  (A). 

The  resonator  parameters  needed  include  the  tube 
quality  factor  (Q) ,  the  speed  of  sound  (C) ,  the  resonance 
frequency  (FO),  the  density  of  the  gas  in  the  tube  (RHO), 
the  resonator  cross-sectional  area  (S),  and  the  equivalent 
length  of  the  resonator  (L) .  For  this  model  the  piston  area 
is  set  equal  to  the  resonator  cross-section  (S).  For  the 
tests  reported  here  they  were  nearly  equal  but  could  differ 
sustantially  in  other  designs, 
b.   Output 

Four  output  plots  are  generated  by  the  graph 
program  that  accompanies  the  model:  (1)  voltage  across  the 
driver  vs.  frequency  (since  we  put  unit  current  into  the 
driver,  this  graph  can  also  be  read  as  the  input  electrical 
impedance),  (2)  piston  velocity  vs.  frequency,  (3)  acoustic 
power  delivered  to  the  load  (resonator)  vs.  frequency,  and 
(4)  electroacoustic  efficiency  vs.  frequency.  Each  of  these 
plots  should  produce  peaks  at  both  the  driver  and  the  tube 
resonant  frequencies . 

The  model  uses  the  electrical  equation  for  the 
coupled  equivalent  circuit  to  calculate  voltage  directly 
(holding  the  current  constant  at  one  ampere)  for  the  first 
output  plot.  It  also  directly  calculates  the  force  on  the 
moving  mass,  F  =  BlI .  From  this  force,  the  piston  velocity 
for   the   second  plot  is  found  by  dividing  the  force  by   the 


mechanical  impedence  (Zm-«=H,).  Pressure  is  calculated  as  tube 
impedence  times  volume  velocity,  and  acoustic  power  is 
pressure  times  volume  velocity  times  the  cosine  of  their 
phase  difference  (Equation  1-8)  for  the  third  plot.  For  the 
final  plot  the  electrical  power  is  calculated  as  voltage 
times  current  times  the  cosine  of  their  phase  difference 
(Equation  1-7),  and  the  model  takes  a  ratio  of  the  electri- 
cal to  the  acoustical  powers  for  the  electroacoustic  effi- 
ciency. Since  the  components  of  these  powers  are  complex  and 
have  phase  differences,  the  actual  powers  are  calculated  as 
one  half  the  real  part  of  [current  times  complex  conjugate 
of  voltage],  and  one  half  the  real  part  of  [pressure  times 
the  complex  conjugate  of  volume  velocity] .  Each  of  the  above 
mentioned  plots  is  the  discussed  output  parameter  versus 
frequency. 

2 .   Assumed  Typical  Parameters 

Table  II-l  lists  the  input  parameters  with  sample 
values  for  the  Altec  290-16K  driver  in  air  with  a  long 
straight  resonant  tube  attached  (see  Chapter  III).  Sample 
outputs  of  the  four  plots  using  these  input  values  are  given 
in  Figures  II-5  through  II-8. 

Figure  II-5  plots  the  driver  voltage  as  a  function 
of  frequency.  Since  the  current  is  held  constant  at  one 
ampere,  this  plot  actually  gives  us  information  on  driver 
electrical  input  impedance  vs.  frequency.  Notice  that  there 
are   two   peaks  in  this  curve.   The  first  peak  occurs   at   a 
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TABLE  II-l.   MODEL  INPUT  VALUES — EXAMPLE  DATA 

Symbol  Meaning  Units  Value 

RM  Mechanical  resistance 

M  Moving  mass 

SSUP  Spring  constant 

RHO  Density  of  gas 

C  Sound  speed 

S  Resonator  X-sec .  area 

BL  Driver  transduction  factor   N/A 

Q  Resonator  quality  factor 

FO  Resonator  resonance  f req .    Hz  240.0 

A  Back  volume  X-sec.  area      m=:  7  qovio 

VOL  Back  volume  m^  1.0X10"- 

LE  Voice  coil  inductance        H  1.0X10-- 

RE  Voice  coil  DC  resistance     H  10.6 

I  Driver  current  A  1.0 

L  Resonator  length  m  0.712 


kg/s 

1  .  33 

kg 

0.0151 

N/m 

9.2X103 

kg/m=' 

1.2 

m/s 

342.0 

m^ 

9.35X10 

N/A 

21  .5 

60  .  0 
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Computer  Model  Sample  Output--Driver 
Voltage  vs.  Frequency  for  One  Ampere 
Current  into  Driver. 
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Figure  II-6. 


Computer  Model  Sample  Output--Piston 
Velocity  vs.  Frequency  for  One  Ampere 
into  Driver. 
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Figure  II-7 


Computer  Model  Sample  Output--Acoustic 
Power  vs.  Frequency  for  One  Ampere. 
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Figure  II-8 


Computer  Model  Sample  Output-- 
Electroacoustic  Efficiency  vs 
Frequency . 
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lower  frequency  and  has  a  much  higher  amplitude  than  the 
second.  This  higher  amplitude  peak  occurs  near  the  driver 
resonant  frequency  ( f  ^  =  { l/27r)  [k/m^]  ^''=2  =  124  Hz).  It 
dominates  the  curve  because  it  is  the  driver  impedance  that 
we  are  plotting,  and  the  voice  coil  motion,  and  consequently 
the  "back  emf , "  is  the  largest  at  that  frequency.  The  fact 
that  there  is  a  second  peak  near  the  tube  resonance  (FO  = 
240  Hz)  shows  the  effect  of  the  coupling  of  the  two  parts  of 
the  system. 

In  the  plot  of  piston  velocity  vs.  frequency  (Figure 
II-6),  the  dominant  peak  again  occurs  near  the  driver 
resonance  for  the  reasons  presented  above.  The  smaller 
coupled  peak  at  the  tube  resonance  is  due  to  the  change  in 
the  load  impedance  which  occurs  near  the  tube  resonance. 

Figure  II-7  shows  the  plot  of  the  acoustic  power  vs. 
frequency.  Now  the  dominant  peak  occurs  at  the  tube 
resonance  because  acoustic  power  is  primarily  a  function  of 
the  tube  parameters .  Here  the  secondary  peak  at  the  driver 
resonance  is  a  consequence  of  the  greatly  increased  piston 
velocity  at  the  driver  resonance  putting  energy  into  the 
smaller  non-resonant  acoustic  load  of  the  tube. 

In  the  plot  of  electroacoustic  efficiency  vs. 
frequency,  Figure  II-8,  the  primary  peak  at  the  tube 
resonance  is  sharp  and  the  secondary  peak  at  the  driver 
resonance  is  small  but  visible.  It  is  important  to  note  that 
the   electroacoustic   efficiency  peaks  at  a  value  of   nearly 


60%  even  though  the  two  individual  resonances  are  separated 
by   nearly  an  octave  in  frequency.   The  peak  acoustic   power 

for  one  amp  as  shown  in  Figure  II-7  is  10.3  W.  The  losses 
(40%)  are  due  primarily  to  Joule  heating  in  the  voice  coil 
DC  resistance  (P^  =  1/2  (I^Rdc)  =  5.3  W) .  This  accounts  for 
34%.  The  remaining  loss  can  probably  be  attributed  to 
mechanical  resistance  losses  in  the  suspension. 

Appendix  A  contains  a  printout  of  the  calculation 
program.  Appendix  B  contains  a  printout  of  the  graph 
program . 
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III.   DRIVER  PARAMETER  MEASUREMENTS 

A.   INTRODUCTION 

In  Sections  II.B  and  II.C  we  discussed  the  choice  of  the 
electrodynamic  transducer  to  drive  the  thermoacoustic  refri- 
gerator. Given  this  choice,  we  needed  to  develop  a  system 
for  measuring  the  mechanical  parameters  of  this  type  of 
driver  to  be  used  in  determining  the  electroacoustic  effi- 
ciency. The  relevant  parameters  (as  discussed  in  Section 
II.D)  are:  (1)  moving  mass,  (2)  suspension  stiffness,  (3) 
transduction  coefficient  (Bl),  and  (4)  mechanical  resis- 
tance. I 

The  driver  we  selected  for  the  preliminary  measurements 
is  the  Altec  model  290-16K  moving-coil  loudspeaker.  There 
are  several  pieces  of  equipment  needed  to  interface  the 
driver  with  the  resonator.  Three  of  these  parts--the  reducer 
cone,  stiffener,  and  dynamic  bellows--become  an  integral 
part  of  the  driver,  and  their  effects  on  the  driver  parame- 
ters must  be  taken  into  account.  We  first  measured  the 
parameters  of  the  driver  alone,  then  added  the  other  parts 
and  determined  the  new  parameters  for  the  entire  system.  The 
driver-bellows  combination  was  then  fitted  to  one  end  of  a 
resonant  tube.  The  actual  resonator  used  in  the  thermo- 
acoustic refrigerator  will  be  a  complex  structure  that  is 
not  easily  modeled.  For  the  purposes  of  the  driver  parameter 
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measurements  we  used  a  resonant  tube  of  uniform  cross- 
section  to  allow  easy  comparison  of  experimental  values  to 
theory.  The  non-driver  end  of  the  closed  tube  housed  several 
microphones.  The  electroacoustic  efficiency  of  the  system  at 
resonance  was  determined  from  measurements  of  electric  power 
delivered  to  the  driver  and  the  acoustic  pressure  at  the 
microphone.  The  acoustic  power  dissipation  was  determined 
from  the  quality  factor  of  the  tube  at  resonance  and  the 
microphone  sensitivity  (the  relationship  between  the  pres- 
sure that  the  microphone  senses  and  its  voltage  output). 

B.   IMPEDANCE  ANALYZER 

The  primary  instrument  used  in  the  testing  and  evalua- 
tion of  the  electrodynamic  driver  was  the  Hewlett  Packard 
Model  •4192A  LF  Impedance  Analyzer.  According  to  the  manufac- 
turers Operation  and  Service  Manual,  the  4192A  is  a  fully 
automatic,  high  performance  test  instrument  designed  to 
measure  a  wide  range  of  impedance  parameters  as  well  as 
gain,  phase,  and  group  delay.  Its  two  measurement  display 
sections  provide  direct  readout  of  the  selected  measurement 
parameters  with  4.5  digit  resolution  along  with  the  appro- 
priate units.  In  normal  mode  operation,  the  4192A  performs 
approximately  five  measurements  per  second.  It  also  provides 
an  average  measurement  mode  (one  measurement  per  second)  for 
higher  resolution  and  noise  rejection.  A  high  speed  mode  can 
perform  ten  measurements  per  second. 
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The  4192A  can  provide  measuring  frequency,  OSC  voltage 
level,  and  DC  bias  voltage  (impedance  measurements  only) 
equivalent  to  actual  operating  conditions .  The  sweep  capa- 
bility of  the  built-in  frequency  synthesizer  and  DC  bias 
source  permits  quick,  convenient,  and  accurate  measurements. 
Measuring  frequency  or  DC  bias  voltage  can  be  automatically 
or  manually  swept  in  either  direction.  OSC  level  can  be 
manually  swept  in  either  direction  in  ImV  increments.  Actual 
test  voltage  across  or  test  signal  current  through  the 
device  under  test  is  also  measured. 

In  amplitude/phase  measurements,  the  4192A  can  measure 
four  transmission  parameters  -  gain/loss,  level,  phase,  and 
group  delay.  In  impedance  measurements,  the  4192A  can 
measure  eleven  impedance  parameters  -  absolute  value  of 
impedance,  absolute  value  of  admittance,  phase  angle,  resis- 
tance, reactance,  conductance,  susceptance,  inductance, 
capacitance,  dissipation  factor,  and  quality  factor.  All 
impedence  measurements  have  a  basic  accuracy  of  0.1%  and  a 
resolution  of  4.5  digits.  Also,  the  circuitry  of  the  4192A 
provides  direct  and  accurate  impedance  measurements  of  both 
grounded  and  floated  devices . 

C.   MECHANICAL  PARAMETERS 

The  most  important  parameters  for  a  moving  coil  trans- 
ducer are  the  moving  mass,  suspension  stiffness,  and  the  Bl 
(transduction)   coefficient.   The  moving  mass  of  the   driver 
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includes  the  mass  of  the  voice  coil,  the  diaphragm,  and  the 
hydrodynamic  mass  of  the  air  being  pushed  by  the  diaphragm. 
The  masses  of  the  reducer  and  stiff ener  are  also  important, 
as  are  the  moving  mass  and  suspension  stiffness  of  the 
bellows.  The  moving  mass  and  stiffness  of  the  combined 
driver-reducer-bellows  system  is  also  important.  The  values 
for  the  combined  system  should  equal  the  sum  of  the  values 
of  the  individual  parts.  The  stiffness,  mass,  and  Bl  coeffi- 
cient give  us  the  information  we  need  to  construct  the 
equivalent  electrical  circuit,  as  discussed  in  Subsection 
II. C. 2.  Also,  we  have  to  insure  that  the  stiffness  of  the 
driver  system  is  just  large  enough  to  produce  a  mechanical 
resonance  at  the  tube  resonant  frequency,  at  which  point  the 
driver  piston  looks  like  a  soft  (or  open)  tube  end.  This  may 
appear  to  be  contradictory  to  having  a  velocity  node  of  the 
standing  wave  at  the  piston,  but  it  is  not.  It  is  how  we  get 
high  efficiency. 

The  important  design  parameters  of  the  tube  are  its 
resonance  frequency  and  quality  factor.  The  length  of  the 
tube  and  the  speed  of  sound  in  the  contained  gas  determine 
the  resonance  frequency. 

1 .   Dynamic  Mass  Loading  of  Driver 

We  used  the  driver's  similarity  to  a  simple  harmonic 
oscillator  to  find  the  moving  mass  and  suspension  stiffness 
of  the  driver  by  taking  advantage  of  the  fact  that  the 
resonance  frequency  squared  equals  the  stiffness  divided   by 
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the  mass.   If  we  add  mass  to  the  driver.   it  should  linearly 
change  the  inverse  of  the  frequency  (or  period)  squared: 

T=  =  l/f=  =  4Tt=(mc  +  m)/k,  (III-l) 

where  m^  is  the  moving  mass,   m  is  the  added  mass,   and  k  is 

the   suspension  stiffness.   The  instrument  we  used  for   this 

measurement    was   the   HP   4192A   Low   Frequency   Impedance 

Analyzer,  described  in  section  B.  We  added  various  masses  of 

putty   (one   to  nine  grams)  to  the  diaphragm  and   found   the 

resonance  frequency  of  the  driver  for  each  mass .   The   putty 

mass   was  measured  on  the  Sartorius  model  2403  mass   balance 

(serial   number   1237).   This  balance  has  an  accuracy   of   ± 

0.00005  g.   However,   the  putty  could  be  expected  to  gain  or 

lose   more   mass  than  that  due  to  handling.   Based   on   mass 

measurements    of   the   putty   after   use,    we   expect    an 

experimental   accuracy  of  ±  0.005  g  for  the  masses.   A   more 

significant  uncertainty  involving  the  masses  arises  from  the 

inability  to  place  the  masses  evenly  on  the  voice  coil  or  to 

produce   a   perfectly  rigid  coupling  between  the  masses   and 

the   voice   coil.   It  is  not  possible  to  assign   a   specific 

numerical  percentage  error  to  this  uncertainty,   but  we  must 

take   it   into  account  when  evaluating  the   data   associated 

with  mass  loading.   A  typical  plot  of  the  period  squared  vs. 

added   mass   in  air  is  shown  in  Figure  III-l.   We   used   the 

least   squares   method  to  fit  a  straight   line   through   the 

data.  We  determined  the  suspension  stiffness  (k  =  •^:?t^/siope) 

to  be: 
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Figure    III-l 


Least  Squares  Fit  of  Period  Squared  vs  . 
Added  Mass  to  Determine  Stiffness  and 
Moving  Mass  of  the  Altec  290-16K  Driver 
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k  =  2.7  ±  0.1  X  lO'*  N/m,  (III-2) 

and  the  moving  mass  (m^  =  intercept/slope) : 

m^  =  5.8  ±  0.1  g,  (III-3) 

for  this  driver  in  air. 

The  moving  mass  of  the  driver  in  vacuum  is   expected 
to   be   less   than  the  value  in  air  because   in   vacuum   the 
driver  is  not  pushing  a  hydrodynamic  mass.   However,   we  did 
not  take  any  measurements  of  this  driver  in  vacuum. 
2 .   Bellows  Parameters 

The  Servometer  custom-designed  dynamic  bellows 
(shown  in  Figure  III-2  with  test  stand)  is  used  to  provide  a 
flexible  gas   seal   between  the  driver  and  the  resonator. 

The  stiffness  of  the  bellows  was  measured  using  two 
independent  methods  -  static  and  dynamic.  For  the  static 
method  we  used  Hooke's  law;  k  =  -F/x.  The  force  was  deter- 
mined by  pulling  a  partial  vacuum  on  the  bellows  and  multi- 
plying the  vacuum  pressure  (read  on  a  manometer)  by  the 
surface  area  of  the  bellows  (A  =  1.1  ±  o.l  X  10-^  m^)  We 
used  a  G.  L.  Collins  Corporation  model  LMS-119V22-01  Linear 
Variable  Differential  Transformer  (LVDT)  (serial  274674) 
connected  to  an  EG&G  128A  Lock-in  Analyzer  to  determine  the 
displacement  of  the  bellows  under  a  partial  vacuum.  Figure 
III-3  shows  the  LVDT  in  its  calibration  stand.  The  calibra- 
tion curve  for  the  LVDT  is  shown  in  Figure  III-4.  The  least 
squares  fit  slope  of  this  curve  gives  us  the  calibration 
constant  for  the  LVDT: 
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Figure    III-2 


Photograph  of  Servometer  Cus torn-Designed 
Dynamic  Bellows — With  and  Without 
Test  Stand. 
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Figure  III-3. 


Photograph  of  Linear  Variable 
Differential  Transformer  (LVDT) 
in  its  Calibration  Stand. 
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LVDT  CALIBRATION 
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dV/dx  =  -0.378  ±  0.0005  V/in  =  -14.9  ±  0.02  V/m.  (III-4) 
Table  III-l  gives  the  static  test  data  for  the  stiffness 
measurement.  By  averaging  the  values  we  obtained: 

k  =  4.3  ±  0.2  X  10^   N/m.  (III-5) 

TABLE  III-l.   BELLOWS  STIFFNESS--STATIC 
MEASUREMENT  RESULTS 

LVDT  Voltage    Pressure    Pressure    Force  (=P.a)   Stiffness 

( mV )         (cm  Hg)     (XlO^Paj        (N)       (XlO'-^N/m) 
14.34  2.9  3. 87  4.39  4.40 

15.62  3.1  4.13  4.69  4.47 

12.89  2.3  3.07  3.48         4.02 

The  dynamic  measurement  was  a  mass  loading  technique 
similar  to  that  used  for  the  driver.  In  this  case,  we 
measured  the  period  as  a  function  of  added  mass  using  an 
Endevco  model  2225M2  accelerometer  (serial  FA80)  and  a 
Nicolet  model  3091  digital  storage  oscilloscope  to  record 
the  free  decay  of  the  bellows  after  it  had  been  disturbed  by 
the  blunt  wooden  end  of  a  cotton  swab.  We  then  plotted 
period  squared  versus  added  mass  and  found  the  stiffness  and 
moving  mass  from  the  slope  and  intercept  of  the  least 
squares  fitted  line  (see  Figure  III-5)  as  discussed  in  the 
previous  Section.  This  gave  us  values  of: 
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k  =  4.56  ±  0.05  X  10^  N/m  (III-6) 

and   m^  =  5.84  ±  0. 05  g.  (III-7) 

This  value  for  moving  mass  includes  the  accelerometer  cable 
and  is  not  representative  of  the  actual  moving  mass  of  the 
bellows.  Subtracting  the  weight  of  the  accelerometer  gives 
us  a  value  of  approximately: 

m^  =  1.5  ±  0.1  g.  {III-8) 

for  the  bellows.  We  know  the  moving  mass  has  to  be  less  than 
2.3  g,  since  that  is  the  bellows'  total  mass--moving  and 
static .  The  values  obtained  for  bellows  stiffness  using  the 
two  different  methods  were  within  experimental  error.  We 
therefore  felt  confident  using  a  value  of  k^  =  4.5  ±  0.1  X 
10^  N/m  in  our  later  calculations. 
3.   Combined  System 

After  we  glued  the  reducer  (m  =  4.8  g),  stiff ener 
plate  {m  =  4.4  g)  and  bellows  to  the  driver,  we  wanted  to 
find  the  stiffness  and  moving  mass  of  the  combined  driver- 
reducer-bellows  system.  We  first  tried  a  dynamic  mass 
loading  measurement  (same  procedure  as  with  the  driver 
alone) ,  but  found  that  the  period  squared  versus  added  mass 
plot  (see  Figure  III-6)  gave  unrealistic  results.  So  instead 
we  used  the  relationship: 

F  =  -kx  =  BlI  (111-9) 

or      k  =  -BlI/x  =  Bl  (di/dV) (dV/dx)     (III-IO) 

and   made  a  static  measurement  using  the  LVDT .   We   expected 

the  additional  equipment  (reducer,   stiffener,   and  bellows) 
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CURRENT  IN  VOICE  COIL  (mA  dc) 

Driver  Displacement  vs.  Voice  Coil 
Current  as  Measured  with  LVDT/ 
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to  change  the  moving  mass  and  suspension  stiffness  of  the 
system  but  not  the  transduction  coefficient  ( Bl ) .  We  applied 
a  DC  current  to  the  driver  and  recorded  the  LVDT  voltage 
output.  then  changed  the  current  and  recorded  the  new  vol- 
tage. This  process  produced  the  values:  ^i  =  283  ±  0.5  mA 
and  AV  =  2.58  ±  0.005  mV .  We  substitute  these  values  and  the 
value  of  Bl  =  21.5  ±  0.5  N/A  (experimental  measurement  of 
this  quantity  is  discussed  in  the  next  section)  into  Equa- 
tion III-IO  for  a  combined  stiffness  value  of: 

k  =  3.5  ±  0.1  X  10^  N/m.  (III-ll) 

If  we  add  the  stiffness  values  for  the  driver  and  bellows 
together  a  total  stiffness  of  k  =  3.2  +  q  ^  y^  j_q-».  ^^^--^  p,-,^. 
the  system  is  obtained.  Thus  the  total  experimental  stiff- 
ness for  the  combined  system  differs  by  about  10%  from  the 
sum  of  the  individual  stiffnesses.  Using  the  relationship 
mo  =  k/(27tf)^,  and  a  measured  resonance  frequency  of 
235  ±  0.5  Hz,  the  combined  moving  mass  is: 

nio  =  16.2  ±  0.5  g,  (III-12) 

compared   to   m^   =  15.9  ±  0.5  g  for  the  individual   masses 

added  together.  The  two  types  of  mass  measurement  agree 
within  2%. 

D.   TRANSDUCTION  COEFFICIENT  ( BL) 

We  determined  the  Bl  coefficient  of  the  driver  using 
three  independent  methods.  The  first  method  was  to  remove 
the   voice   coil  from  the  magnet  and   directly   measure   the 
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magnetic  induction  using  a  Dowdy  RFL  Industries  model  912 
Guassmeter  (serial  808)  using  a  Hall  effect  probe  model 
912015  (serial  21122)  (B  =  2.2  ±  0.05  Tesla).  We  then 
measured  1  (the  length  of  the  wire)  by  counting  the  number 
of  turns  of  wire  in  the  coil  (under  a  microscope)  and 
multiplying  this  number  (4-4)  by  the  circumference  of  the 
coil  (1  =  44  X  7t  X  diameter  =  9.98  ±  0.005  m).  Combining 
these  two  measurements  gave  us  the  value: 

Bl  =  22.0  ±  0.7  N/A.  (III-13) 

In  the  second  method  we  used  the  relationship: 
F  =  BlI  =  -kx  or  Bl  =  -kx/I  =  -k (dx/dV) (dV/dl ) .  (III-14) 
We  previously  determined  k  for  the  driver  alone  (Equation 
III-2),  and  now  used  the  LVDT  and  lock-in  analyzer  to  find 
dx/di.  The  LVDT  calibration  was  discussed  in  the  previous 
Section,  with  a  calculated  value  of  dV/dx  =  -14.9±o.05  V/m. 
We  then  measured  the  diaphragm  displacement  as  a  function  of 
DC-current  through  the  voice  coil.  This  gave  us  a  value  of 
dV/dl  =  11.75  ±  0.05  mV/A  (see  Figure  III-7).  Using  the 
relation  in  Equation  III-14  we  obtain  the  value: 

Bl  =  21.2  ±  0.1  N/A.  (III-15) 

In  the  third  method  we  used  the  relationship: 

F  =  Bli  =  mg  or  Bl  =  mg/i.  (III-16) 
For  this  we  measured  the  change  in  DC-current  through  the 
coil  as  a  function  of  added  mass,  using  the  LVDT  as  a 
monitor  to  keep  the  displacement  of  the  diaphragm   constant. 
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Thi"se   diflerent  mass /current:  combinations  { s  •=:<=:  Tdblt:  1 1 1-2) 
produced  an  average  value  of: 

Bl  =  21.4  ±  0.5  N/A.  (III-17) 

Averaging  the  results  of  these  three  methods  gave  us  a 
value  of  Bl  =  21.5  ±  0.5  N/A  with  2%  deviation  among  the 
three  methods . 

TABLE  III-2.   STATIC  MEASUREMENT  OF  ALTEC  290-16K 
TRANSDUCTION  COEFFICIENT  (BL) 


Am  ( g ) 
132  .  5 
157.5 
257.5 


111     (mA) 
64.4 
74.0 
108.2 


Bl  =  mg/I  (N/A) 
20.2 
20  .  9 
23.3 


E.   MICROPHONE  CALIBRATION 

Microphone  calibration  is  an  important  part  of  the 
process  leading  to  the  measurement  of  the  electroacoustic 
efficiency.  A  microphone  is  used  to  measure  the  acoustic 
pressure  inside  the  resonator  by  converting  the  acoustic 
pressure  to  an  output  voltage.  The  ratio  of  the  open  circuit 
voltage  to  the  pressure  is  called  the  open  circuit 
microphone  sensitivity. 

We  used  three  microphones  for  the  calibration:  (1)  a 
General    Radio   (GR)   model   1962-9611   half-inch   electret 
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condenser  microphone  (serial  11811)  with  GR  preamplifier 
model  1560-P42  (serial  4061);  (2)  an  Endevco  model  8510B-5 
piezoresistive  bridge  microphone  (serial  73MJ);  and  (3)  an 
1/8  inch  B&K  model  4138  microphone  (serial  19683) .  We  first 
calibrated  the  GR  using  a  GR  type  1562-A  sound-level  cali- 
brator (serial  19683,  114dB  SPL ,  250  Hz),  and  obtained  a  GR 
sensitivity  of  M^,  =  7.42  mV/Pa.  This  compares  within  about 
l7o  of  the  calibration  certificate  provided  by  the  manufac- 
turer of  7.50  mV/Pa.  Additional  checks  were  made  with  a  B&K 
pistonphone  calibrator  which  substantiated  the  GR  calibra- 
tion results. 

To  calibrate  the  other  two  microphones  we  attached  the 
resonator  tube  to  the  290-16K  driver.  We  used  an  HP  3325A 
Synthesizer/Function  Generator  to  excite  the  driver  at  240 
Hz,  the  resonant  frequency  of  the  tube.  We  recorded  the 
output  voltages  for  all  three  microphones  and  used  the 
comparison  relationship: 

Mi/Vi  =  M^/\j^  (III-18) 

to  determine  the  sensitivities  of  the  B&K  and  Endevco.  The 
B&K  sensitivity  was  determined  to  be  Mb  =  0.463  mV/Pa, 
comparing  within  2%  of  the  manufacturer's  specification 
value  of  0.470  mV/Pa  (adjusted  to  include  the  capacitance  of 
the  B&K  model  VA0160  1/8"  to  1/4"  adapter  with  an  expected  4 
dB  "line  loss").  The  sensitivity  of  the  Endevco  microphone 
was  calculated  to  be  8.17  X  10-*  V/Pa ,  comparing  within  1% 
of  the  manufacturer's  specification  of  8.09  X  10-«-  y/Pa   For 
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all  three  microphones  we  measured  sensitivity  values  within 
2%  of  the  manufacturers'  specifications. 

F.   ELECTROACOUSTIC  EFFICIENCY 

The  electroacoustic  efficiency  is  the  ratio  of  the 
acoustic  power  dissipated  in  the  resonator  to  the  electric 
power  delivered  to  the  driver.  The  electric  power  was 
determined  using  the  relationship: 

P-x,<=  =  V  I  cos  e,  (III-19) 

where  V  is  the  voltage  across  the  driver  terminals,  I  is  the 
current  flowing  in  the  driver,  and  e  is  the  phase  angle 
between  the  voltage  and  the  current.  The  values  for  voltage, 
current  and  phase  angle  delivered  to  the  driver  were 
obtained  using  the  HP  4192A. 

Calculation  of  the  acoustic  power  was  a  little  more 
complex.  We  started  with  the  basic  relationship  for  the 
quality  factor  of  the  tube: 

Q  =  27tE.3-r/E,_c ,  (III-20) 

where  Eq-t  is  the  energy  stored  in  the  tube  and  Ei_c  is  the 
energy  lost  per  cycle.  The  acoustic  power  generated  in  the 
tube  is  the  energy  lost  per  cycle  times  the  number  of  cycles 
per  second,  or  frequency: 

P«c  =  ELcf  =  27tfEsx/Q.  (III-21) 

Following   Burmaster   (1985)  the  energy  stored   in   a   plane 

standing    sound   wave   is   equal   to   the   energy    density 

integrated  over  the  volume  of  gas  in  the  tube  (VOL): 
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Est  =  (P,-m.^/2;5P.=  )  »VOL,  (III-22) 

where  P^„.  is  the  root  mean  square  acoustic  pressure  in  the 
tube,  V  ^2  the  ratio  of  the  specific  heats  of  the  gas  in  the 

tube  (1.4  for  air),  and  P^  is  the  ambient  pressure  (one 
atmosphere  or  1.013  X  IC  Pa  for  these  measurements).  The 
r.m.s.  acoustic  pressure  is  the  r.m.s.  microphone  output 
voltage  divided  its  sensitivity  (V/M)  as  discussed  in  the 
previous  Section.  The  expression  for  the  acoustic  power 
becomes : 

F«c:  =  (JT-fV=VOL)/(QM  =  -iP^)  .  (III-23) 

The   only   remaining  quantity  in   Equation   III-23   that 

needs  to  be  calculated  separately  from  the  efficiency 
measurement  is  the  quality  factor  of  the  tube  (Q) .  To 
determine  Q  we  removed  the  290-16K  driver  from  the  resonator 
and  replaced  it  with  an  electret  transducer  to  drive  the 
tube.  This  was  necessary  because  there  were  losses  due  to 
the  290-16K  that  made  the  tube  Q  appear  to  be  lower  than  it 
actually  is  (we  discovered  this  problem  while  investgating 
the  reasons  for  the  calculated  efficiency  to  be  greater  than 
one).  We  used  the  B&K  microphone  for  the  pressure  (converted 
to  voltage)  readings  at  the  other  end  of  the  tube.  We  plot- 
ted voltage  as  a  function  of  frequency,  and  determined  the 
resonant  frequency  (f^)  and  frequencies  of  the  half  power 
points  (f^  and  fi).  Using  the  relationship: 

Q  =  f^/(f^-f^)  (111-2^) 

we  obtained  a  value  of  Q  =  60  ±  l  for  the  tube. 
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To  determine  the  electroacoustic  efficiency  we  returned 
the  tube  to  the  290-16K  driver  (see  Figure  III-8)  and  drove 
the  tube  at  resonance.  We  recorded  the  voltage,  current,  and 
phase  angle  for  the  driver  and  the  B&K  output  voltage.  The 
electroacoustic  efficiency  was  then  calculated  using: 

nE«  =  P«c/P-x-=.  (III-25) 

Table  III-3  lists  the  electric  power,  acoustic  power,  and 
efficiency  for  five  test  runs  at  driver  currents  ranging 
from  3-38  mA  (with  corresponding  voltages  0.15-1.43  V).  The 
electroacoustic  efficiency  ranged  from  BA . 2  for  the  lowest 
driver  current  to  75.9  for  the  highest. 

TABLE  III-3.   ALTEC  290-16K  ELECTROACOUSTIC 

EFFICIENCY  MEASUREMENTS 

He^  ( % ) 
84.2 

79.2 

77.9 

77.2 

75.9 


P-x-<=    (mW) 

P«c 

(mW) 

0.56 

0. 

,47 

3.41 

2. 

.70 

10.0 

7. 

,79 

17.0 

13. 

.1 

54.2 

41. 

,2 

G.   COMPARISON  TO  COMPUTER  MODEL 

The  final  task  in  the  evaluation  procedure  for  the 
electrodynamic  driver  is  to  compare  the  measured  electro- 
acoustic  efficiency  to  that  predicted  by  the  computer  model 
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Figure  III-8. 


Photograph  of  Altec  290-16K  Driver 
Assembly  with  Resonator  Tube  and 
Endevco  Microphone. 
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discussed  in  Section  II. D.  We  have  already  measured  or 
calculated  all  the  input  values  for  the  model  except  for  the 
driver's  mechanical  resistance  (R„)_ 

1 .   Mechanical  Resistance 

We  can  determine  Rm  by  using  the  relationship: 

Rm  =  27rf^m<=./Q,  (III-26) 

where  f^  is  the  resonant  frequency,  m^  is  the  moving  mass 
and  Q  is  the  quality  factor  of  the  driver.  There  are  several 
methods  available  for  determining  the  driver's  quality 
factor.  One  is  to  plot  the  magnitude  of  the  modified  driver 
impedance  vs.  frequency  and  pick  the  half  power  points  off 
the  graph.  Then,  as  with  the  tube,  Q  is  given  by: 

Q  =  fc/{fc.-fi  )  ,  (III-27) 

where  f^  and  fi  are  the  frequencies  associated  with  the 
upper  and  lower  half  power  points  (located  at  [Zl^.^/V^", 
where  [Z]  is  the  magnitude  of  the  driver  impedance  which  is 
modified  by  subtracting  the  direct  current  (DC)  resistance 
from  the  measured  complex  irapadance). 

Another   method  is  to  use  the  relationship   [Hofler, 
1986]: 

Q  =  J-  f^  dO/df  ,  (III-2S) 

where  t>  is  the  phase  angle  between  driver  resistance  (R)  and 
reactance  (X),  the  two  components  of  modified  impedance  (Z  = 
R  +  iX,  where  R  is  the  measured  resistance  minus  the  DC 
(direct  current)  resistance).  This  modified  impedance  is 
proportional  to  the  voice  coil  velocity.   The  quantity  dt>/df 
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is  the  slope  of  txhe  phase  angle  versus  frequency  curve 
evaluated  at  the  resonant  frequency  (where  the  curve  is  a 
straight  line) .  We  used  the  HP4192A  to  obtain  the  data  for 
these  measurements . 

Figure   III-9  shows  a  plot  of  modified  impedance  vs. 
frequency  for  the  Altec  290-16K  in  air.   From  this  graph   we 

find  fc  =  235.2  Hz,  fo.=249.6  Hz  and  fi  =222.0  Hz.  The  Q  for 
these  values  is  8.52.  Figure  III-IO  shows  a  plot  of  phase 
angle  vs.  frequency  for  the  same  data  set.  Over  this  large 
range  of  frequencies  the  phase  angle  curve  approximates  an 
arctangent.  At  resonance,  however,  where  dy/df  is  evaluated, 
the  curve  is  linear  (see  Figure  III-ll).  Equation  III-28  is 
the  slope  of  the  arctangent  at  resonance  (fo).  The  slope  of 
the  plotted  data  is  71.8X10-='  sec  (dt>  was  converted  from 
degrees  to  radians).  The  Q  calculated  from  this  method  is 
8.44.  These  values  for  Q  agree  within  1%.  From  these  values 
we  get  R„=2.77  kg/s . 

We  tried  yet  another  method  for  determining  driver 
Q.  This  third  method  measures  the  decay  constant,  x,  using  a 
digital  storage  oscilloscope  to  record  the  "ring  down"  of 
the  driver.  This  technique  takes  less  time  to  run  the  exper- 
iment and  reduce  the  data  using  the  relationship: 

Q  =  a>^/2  or  X  =  2Q/aj  =  QT/TT  (III-29) 
where  a>  is  the  resonant  angular  frequency  in  radians/sec,  T 
is   the   period  in  seconds,   and  t  is  the   free   decay   time 
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Figure    III-9 


Altec  290-16K  Modified  Impedance 
Magnitude  vs.  Frequency. 
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Figure  III-IO.   Altec  290-16K  Modified  Impedance 

Phase  Angle  vs.  Frequency — Full 
Frequency  Range. 
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Least  Squares  Fit  of  Altec  290-16K 
Modified  Impedance  Phase  Angle  as  a 
Function  of  Frequency  Near  Resonance 
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constant   in   seconds.   This   method  will   be   discussed   in 
greater  detail  in  Chapter  IV. 
2 .   Model  Results 

Table  III-4  lists  the  computer  model  input  values 
for  the  Altec  290-16K.  Figures  III-12  through  III-15  show 
the  output  plots  for  these  values  as  discussed  in  Section 
II. D.  In  this  case  the  resonance  peaks  are  much  closer 
together,  and  the  magnitude  of  the  secondary  peak  is  closer 
in  value  to  the  primary  peak.  Once  again  the  driver  reso- 
nance is  at  a  lower  frequency  than  the  resonator  resonance. 
It  is  interesting  to  note  that  the  modeled  efficiency  (about 
607o)  is  15-20%  less  than  the  measured  efficiency  reported  in 
the  previous  Section.  We  believe  that  this  results  from  the 
fact  that  for  the  computer  model ,  a  high  tube  Q  means  that 
the  driver  losses,  represented  by  R„ ,  play  a  more  important 
role  in  the  efficiency  calculation,  and  the  results  are  very 
sensitive  to  small  errors  in  R„ .  The  primary  ways  for  the 
driver  energy  to  be  dissipated  are:  (1)  Joule  heating  due  to 
DC  resistance  in  voice  coil  and  eddy  currents  induced  in  the 
magnet  structure  due  to  the  time-varying  magnetic  fields 
generated  by  the  AC  voice  coil  current  (small  but  not  negli- 
gible) ;  (2)  mechanical  resistance  of  driver  (R,„);  and  (3) 
energy  transferred  to  acoustic  load.  Another  possible  reason 
for  the  different  efficiencies  is  that  the  efficiency  is  not 
independent  of  current  hence  the  real  physical  system  has 
some   non-linearity,   as   Table   III-3   clearly   shows.   The 
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TABLE  III-4.   COMPUTER  MODEL  INPUT--ALTEC  290-16K  DRIVER 

Symbol  Meaning  Units         Altec 

RM  Mechanical  resistance 

M  Moving  mass 

SSUP  Spring  constant 

RHO  Density  of  gas 

C  Sound  speed 

S  Resonator  X-sec.  area 

BL  Driver  trans,  factor 

Q  Resonator  quality  factor 

FO  Resonator  resonance  f req . 

A  Back  volume  X-sec .  area 

VOL  Back  volume 

LE  Voice  coil  inductance 

RE  Voice  coil  DC  resistance 

I  Driver  current 

L  Resonator  length 


kg/s 

2.77 

kg 

0.0162 

N/m 

3.52X10' 

kg/m^ 

1  .2 

m/s 

342  .  0 

m= 

9. 35X10 

N/A 

21  .5 

60.0 

Hz 

240.0 

m=^ 

4 . 13X10 

m"^ 

2. 3X10" 

H 

1 .0X10- 

Pl 

10.6 

A 

1  .0 

m 

0.712 
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Figure  III-12 


Computer  Model  Output  for 
Altec  290-16K--Driver  Voltage 
vs.  Frequency  for  One  Amp. 
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Figure  III-13 


Computer  Model  Output  for 
Altec  290-16K--Piston  Velocity 
vs.  Frequency  for  One  Amp. 


82 


Figure  III-14 


Computer  Model  Output  for 
Altec  290-16K--Acoustic  Power 
vs.  Frequency  for  One  Amp. 
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Figure  III-15. 


Computer  Model  Output  for 
Altec  290-16K--Electroacoustic 
Efficiency  vs.  Frequency. 


84 


computer  model  is  based  on  a  linear  theory  and  vjill  have   an 
efficiency  which  is  independent  of  drive  current. 

H.   DRIVER  RESONANCE  ALTERATION 

The  refrigerator  must  operate  at  the  resonant  frequency 
of  the  resonator  to  obtain  the  highest  electroacoustic  effi- 
ciency and  to  maintain  stack  position  with  respect  to  the 
standing  acoustic  wave  {discussed  in  Sub-section  II. A. 4).  It 
appears  from  the  computer  model  results  and  our  experiments 
so  far  that  the  operating  frequency  should  also  be  near  the 
driver's  resonance  to  get  the  most  power  out  of  the  driver. 
The  resonant  frequency  of  the  resonator  is  determined  by  the 
effective  length  of  the  resonator  and  the  speed  of  sound  in 
the  contained  gas.  The  effective  length  for  a  tube  of  uni- 
form cross-section  and  closed  end  is  the  actual  measured 
length.  For  the  STAR  resonator,  however,  the  effective 
length  is  the  length  that  would  produce  the  measured  reso- 
nant frequency  for  the  given  sound  speed  if  the  resonator 
were  a  uniform  cross-section  closed-ended  tube  (see  Figure 
1-2  for  approximate  shape  of  resonator).  The  sound  speed  can 
be  altered  by  changing  the  percentage  of  xenon  in  the 
helium-xenon  mixture,  since  the  sound  speed  depends  on  the 
density  of  the  mixture,  and  the  density  depends  on  the 
percentages  of  the  individual  gases  present  (see  Equations 
IV-24  and  IV-26).  The  gas  mixture  and  effective  resonator 
length   will   be   chosen   to   maximize   the   coefficient   of 


performance  of  the  refrigerator.  Therefore  if  we  want  to 
match  the  resonances,  we  will  have  to  alter  the  resonance  of 
the  driver  once  the  resonant  frequency  of  the  resonator  has 
been  determined . 

The  resonant  frequency  of  the  driver  is  proportional  to 
the  square  root  of  the  stiffness  divided  by  the  moving  mass. 
If  we  want  to  raise  the  resonant  frequency  of  the  driver  we 
can  either  increase  the  stiffness  or  decrease  the  mass. 
Conversely,  to  lower  the  resonant  frequency  we  can  either 
decrease  the  stiffness  or  increase  the  mass.  Generally  we 
prefer  to  decrease  stiffness  and  mass  to  get  a  larger 
velocity.  Reduction  of  stiffness  usually  also  reduces  the 
mechanical  resistance. 

For  the  Altec  290-16K  attached  to  the  uniform  cross- 
section  resonant  tube  in  these  preliminary  experiments ,  the 
frequencies  were  already  well  matched  (note  the  closeness  of 
the  peaks  in  Figures  III-12  through  III-14)  and  changing  the 
driver  frequency  was  not  required.  However,  we  know  that  in 
future  applications  we  may  need  to  alter,  most  probably 
lower,  the  driver  resonant  frequency.  We  developed  an 
experimental  procedure  to  do  this. 

Adding  moving  mass  to  the  driver  would  appear  to  be  the 
easiest  way  to  lower  the  driver  resonance.  Unfortunately  the 
added  mass  would  also  reduce  the  velocity  by  increasing  the 
mecahnical  impedance  above  resonance.  We  therefore  chose  to 
experiment  with  decreasing  the  driver  stiffness.  Recall  that 


86 


the  driver-reducer-bellows  system  stiffness  includes  the 
stiff fness  of  the  bellows  (4.5  X  lO-"^  N/m)  and  the  suspension 
stiffness  of  the  driver  phenolic  surround  (29  X  10"^  N/m). 
The  driver  suspension  was  the  obvious  choice  for 
experimental  alteration  for  two  reasons:  (1)  it  is  over  six 
times  as  large  as  the  bellows  stiffness,  and  (2)  it  would  be 
very  difficult  for  us  to  change  the  bellows  stiffness 
without  affecting  its  performance. 

Since  the  290-15K  voice  coil  surround  is  made  of  a 
phenolic  material,  it  was  relatively  easy  to  cut  through  it 
with  an  ordinary  scalpel.  We  made  a  series  of  successive 
cuts  in  the  surround  material  and  ran  measurements  of  the 
driver  resonant  frequency  after  each  cut.  The  results  of 
this  experiment  are  given  in  Table  III-5.  We  plotted  the 
phase  angle  between  the  driver's  electromotive  force  (emf) 
and  the  voice  coil  displacement  as  a  function  of  frequency. 
For  the  emf  we  measured  the  voltage  across  the  driver  termi- 
nals in  channel  A  of  the  HP  4192A  Impedance  Analyzer  and 
assumed  that  the  emf  was  in  phase  with  the  driver  current 
and  hence  force.  We  used  the  photonic  sensor  to  measure  the 
displacement  by  connecting  its  voltage  output  to  channel  B. 
The  4192A  then  displayed  the  phase  angle  between  the  two 
voltages  representing  force  and  displacement.  The  driver  is 
at  resonance  when  this  phase  angle  equals  90  degrees .  The 
quality  factor  (Q)  is  equal  to  one-half  the  resonant  fre- 
quency  times   the  slope  of  the  phase   angle   vs.   frequency 
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curve  evaluated  at  resonance  (see  Equation  III-28).  Using 
the  previously  determined  values  for  moving  mass  (16  g)  and 
bellows  stiffness,  we  calculated  the  suspension  stiffness: 

^m^-»     =     kfcot     -     kfc>«iiow«  (III-30) 

where  ki^ct  =  '^n^f^^m^.  We  calculated  the  mechanical 
resistance   (Rm)  using  equation  III-26.   We  stopped   cutting 

the   surround  when  its  stiffness  was  approximately  equal   to 

that  of  the  bellows  and  therefore  the  total  stiffness  was  no 

longer  suspension  dominated. 

TABLE  III-5.   RESULTS  OF 
DRIVER  RESONANT  FREQUENCY  ALTERATION 

Percent  Surround   Resonant   Quality   Suspension   Mechanical 
Remaining     Frequency   Factor    Stiffness    Resistance 


(approx 

•) 

(Hz) 

(XIO^ 

^N/m) 

(kg/s) 

100 

241  .  0 

8. 

.55 

32  .  2 

2.83 

100  (six 

3 

lits) 

229.5 

8. 

.09 

28.8 

2.85 

80 

213.5 

5. 

.  17 

24.3 

4.  15 

50 

183.5 

8. 

.65 

16.8 

2.  13 

20 

119.5 

8. 

.56 

4.52 

1  .40 

Despite  the  anomoly  at  the  80%  data  point,  we  feel  that 
this  experiment  shows  that  partial  removal  of  the  surround 
can  reduce  stiffness,  resonant  frequency.  and  mechanical 
resistance.   It  appears  from  the  above  data  that  the  flexing 
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of   the  phenolic  surround  is  responsible  for  the  majority  of 
the  mechanical  losses  in  the  driver. 
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IV.   STAR  DRIVER  PARAMETERS 

A.   INTRODUCTION 

This  chapter  discusses  the  selection  of  the  electro- 
dynamic  driver,  describes  the  parameter  measurements  on  the 
chosen  driver,  and  estimates  the  electroacoustic  efficiency 
of  the  STAR  using  the  computer  model  described  in  Section 
II. D. 

1 .   Interface  with  NASA  GAS  Cannister 

In  addition  to  the  maximum  weight  of  200  pounds 
discussed  in  Section  I.B,  restrictions  due  to  the  Get  Away 
Special  (GAS)  cannister  dimensions  give  us  a  maximum  payload 
height  (28.25  inches)  and  diameter  (19.75  inches).  In  our 
alloweu  height  we  must  fit  two  battery  layers,  the  driver 
with  its  housing,  the  resonator  assembly,  and  its  associated 
vacuum  space  and  thermal  insulation.  The  battery  layers  are 
each  about  four  inches  high,  for  a  total  battery  height  of 
eight  inches.  This  leaves  about  20  inches  for  the  driver  and 
resonator  assemblies  (see  Figure  1-4  for  approximate 
configuration) .  The  electronics  package  and  gas  reservoir 
should  fit  in  the  area  radially  surrounding  the  resonator, 
and  will  therefore  not  add  to  the  height.  Although  the 
diameter  of  the  cannister  volume  is  19.75  inches,  the  driver 
housing  must  attach  to  the  standard  GAS  can  lid  at  a  9.5 
inch   diameter   bolt  circle,   otherwise  we  would  need  a   lid 
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specially  designed  for  our  project.  As  far  as  weight  is 
concerned.  approximately  every  pound  we  save  elsewhere 
allows  us  to  carry  an  extra  battery  cell,  which  results  in 
an  extra  ten  watt-hours  of  available  energy. 

2 .  Maximizing  Electroacoustic  Efficiency 

In  addition  to  the  space  and  weight  restrictions,  we 
want  to  maximize  the  electroacoustic  efficiency  of  the 
driver.  We  can  do  this  by  obtaining  the  maximum  value  for 
the  combination  of  parameters  Bl/R„mo.  A  larger  Bl  will 
give  us  a  larger  driving  force  (F=Bli),  which  in  turn  gives 
us  a  larger  acoustic  volume  velocity,  one  of  the  components 
of  acoustic  power.  For  a  given  value  of  magnetic  induction 
(B),  a  longer  voice  coil  (achieved  by  having  a  larger  voice 
coil  diameter  with  the  same  number  of  turns,  or  the  same 
diameter  and  more  turns)  will  result  in  a  larger  Bl .  This 
longer  1  would  also  result  in  a  greater  voice  coil 
electrical  resistance,  so  the  two  effects  cancel.  We  also 
want  to  minimize  the  moving  mass  (m^)  to  allow  greater 
velocities  and  minimize  the  driver's  mechanical  resistance 
(Rm)  in  order  to  minimize  mechanical  losses  in  the  driver. 

3 .  Driver  Selection 

We  had  two  off-the-shelf  drivers  available  for  con- 
sideration. The  Altec  260-16K  was  discussed  in  Chapter  III. 
The  other  driver  was  a  JBL  2445J .  Both  drivers  were  designed 
for  applications  somewhat  different  from  the  STAR.  Figure 
IV-1   shows   a  side  view  of  the  2'4A5J   driver.   The   hatched 
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Figure  IV-1 


Schematic  Drawing  of  JBL  2'^A5J  Driver 
Hatched  Area  Shows  Unnecessary  Throat 
Section. 
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area,  known  as  the  throat  of  the  driver,  is  considered 
unnecessary  for  our  application,  and  adds  a  lot  of  extra 
weight  to  the  driver.  The  Altec  290-16K  has  a  similar  excess 
throat  volume.  We  contacted  JBL  and  asked  if  they  would 
custom-design  a  driver  which  would  eliminate  this  excess 
weight.  Fortunately  for  our  project,  JBL  was  in  the  process 
of  designing  a  new  driver,  the  2450 J,  which  uses  neodymium- 
iron-boron  magnets  instead  of  the  traditional  ceramic.  They 
agreed  to  custom-design  a  version  of  the  2450 J  for  our  use. 
The  resulting  driver  (hereafter  referred  to  as  NIB)  is  shown 
in  Figure  IV-2.  The  NIB  driver  accepts  the  same  voice  coils 
(four  inch  diameter,  aluminum  wire,  and  titanium  suspen- 
sion/diaphragm) as  the  2445J  (ceramic  magnet)  driver.  JBL 
claims  that  it  should  also  have  the  same  properties, 
including  Bl . 

In  the  remainder  of  this  chapter  we  describe  the 
parameter  tests  on  the  NIB  driver.  We  then  compare  the  NIB 
driver  to  the  290-16K  driver.  The  results  of  the  2445J 
parameter  measurements  are  given  for  comparison.  Finally  we 
present  results  generated  by  the  computer  model  for  the  NIB 
driver  given  various  acoustic  load  impedances. 

B.   NIB  DRIVER  PARAMETERS 

For  the  measurement  of  the  NIB  driver  parameters,  we 
used  the  techniques  we  developed  in  Chapter  III  for  the 
measurement  of  the  290-16K  driver  parameters. 


Figure  IV-2 


Photographs  of  JBL  Neodymium- 
Iron-Boron  Driver. 
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1 .   Moving  Mass  and  Suspension  Stiffness 

The  moving  mass  and  suspension  stiffness  were 
calculated  using  the  mass  loading  technique.  We  placed  small 
masses  (0.5  to  3.5  grams)  of  black  putty  on  top  of  the 
voice  coil  in  half-gram  increments.  For  each  added  mass  the 
resonant  frequency  of  the  driver  was  found  using  the  HP 
4192A  LP  Impedance  Analyzer  (see  Section  III.B  for  a 
description  of  this  instrument).  Figure  IV-3  shows  a  sample 
plot  of  period  squared  vs.  added  mass  for  the  NIB  driver  in 
air.  The  stiffness  was  found  by  calculating  four  pi  squared 
divided  by  the  slope.  The  moving  mass  is  equal  to  the 
intercept  divided  by  the  slope.  For  the  NIB  in  air  these 
values  are: 

k    =  9.8  ±  0.1  X  10*  kg/s^  (IV-1) 

and      m^  =  7.8  ±  0.1  g,  (IV-2) 

with   a   resonant   frequency  of   564   Hz.   We   repeated   the 
measurement  in  a  vacuum  chamber  and  found: 

k  =  8.6±o.lX  10-»  kg/s==  (IV-3) 

and      m^  =  3.2  ±  0.1  g,  (IV-4) 

with  a  resonant  frequency  of  833  Hz  (see  Figure  IV-4). 

The  moving  mass  is  greater  in  air  because  of  the 
driver  having  to  move  the  air  itself.  The  discrepancy  in  the 
stiffness  is  not  as  easy  to  explain.  To  determine  the  actual 
value  of  k,  we  used  the  static  mass  loading  method.  This 
gave  us  a  value  of: 
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NIB  MASS  LOADING  -  AIR 
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k  =  6.4  ±  0.2  X  10-*  kg/s=.  (IV-5) 

(The   static   loading  measurements  will  be  discussed  in   the 

next  Sub-section.) 

This  last  value  for  suspension  stiffness  differs  by 
about  40%  from  the  first  two.  We  did  not  have  sufficient 
time  to  test  this  driver  until  we  were  sure  we  had  the 
"correct"  value  for  the  stiffness.  We  recommend  that  the 
follow-on  students  do  so.  Lacking  the  correct  value,  we  ran 
the  computer  model  for  two  sets  of  data:  one  set  represen- 
ting the  "high"  k  (9.0  X  10^  kg/s==)  and  one  using  the  "low" 
k  (6.4  X  10-*  kg/s=2)  . 

2  .   Transduction  Coefficient  ( Bl ) 

The  transduction  coefficient  (Bl)  was  calculated 
using  two  methods.  The  first  method  was  to  measure  B  (magne- 
tic induction  in  the  voice  coil  gap)  and  1  (length  of  the 
voice  coil)  directly.  We  measured  B  with  the  Dowty  RFL 
Industries  model  912  Guassmeter  (serial  303)  using  a  Hall 
effect  probe  model  912015  (serial  21122) .  The  maximum 
reading  using  this  probe  was  18.7  kGauss  =  1.87  Tesla  (T). 
This  value  differs  by  less  than  one  percent  from  the  value 
stamped  on  the  driver  by  the  manufacturer:  1.858  T.  We  will 
therfore  use  a  value  of  B=1.86  ±  0.01  T  for  this  driver.  We 
then  counted  the  turns  of  aluminum  wire  on  the  voice  coil  (n 
=  30  ±  1  turns)  and  calculated  1  using  the  relationship: 

1  =  7f  d  n.  (IV-6) 
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For  the  9.98  -  0.005  cm  diameter  coil  we  calculated  a  value 
of  1  =  9.4  ±  0.3m.  Multiplying  these  two  values  together 
gives  us  for  the  transduction  coefficient: 

Bl  =  17.5  ±  0.6  N/A.  {IV-7) 

The   other   method   used  to  determine   Bl   used   the 
relationship : 

F  =  BlI  =  -kx  or  Bl  =  -kx/I ,  {IV-8) 
where  I  is  the  direct  current  applied  to  the  coil,  k  is  the 
spring  constant  measured  earlier,  and  x  is  the  displacement 
of  the  voice  coil  producea  by  the  current.  For  this  measure- 
ment on  the  290-16K  driver  we  had  used  a  Linear  Variable 
Differential  Transformer  (LVDT).  For  the  NIB  driver, 
however,  we  used  Mechanical  Technology  Incorporated ' s  MTI- 
1000  photonic  sensor  (serial  38010818)  with  plug-in 
cartridge  model  MTI-3S12  (serial  1366)  and  probe  model  125R. 
Figure  IV-5  shows  the  photonic  sensor  with  the  NIB  driver. 
One  major  advantage  of  the  photonic  sensor  over  the  LVDT  is 
that  the  photonic  does  not  need  to  make  physical  contact 
with  the  surface  it's  measuring  (it  uses  light)  whereas  the 
LVDT  does.  With  the  Altec  290-16K's  rigid  phenolic  dia- 
phragm, the  contact  with  the  LVDT  probably  did  not  affect 
the  measurement.  The  NIB  diaphragm,  however,  is  made  of  thin 
titanium  and  is  quite  flexible.  It  would  have  been  impos- 
sible to  get  an  accurate  reading  of  the  displacement  of  this 
diaphragm  using  the  LVDT.  We  calibrated  the  photonic  sensor 
using   calibration  fixture  KD-CH-III  A.   The  output   voltage 
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Figure  IV-5.   MTI-1000  Photonic  Sensor  with  NIB  Driver 
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from  the  sensor  was  read  on  a  HP  3473A  multimeter  (serial 
A110136).  We  then  plotted  DC  voltage  vs.  displacement  in 
millimeters.  The  resulting  graph  is  shown  in  Figure  IV-6. 
The  slope  of  this  line  gives  us  the  calibration  constant: 

dV/dx  =  -2.086  ±  0.001  V/mm.  {IV-9) 

This  value  is  within  two  percent  of  the  calibration  value 
provided  by  the  manufacturer  of  dV/dx  =  2.045  V/mm  (back 
slope) . 

The  next  step  was  to  run  DC  current  through  the  coil 
and  record  the  displacement  of  the  diaphragm  (via  a  voltage 
reading  from  the  photonic  sensor).  For  the  voltage  supply  we 
used  the  HP  467A  power  amplifier  (serial  444-00243).  We 
found  the  current  by  measuring  the  voltage  (with  the  HP 
3476A   digital  multimeter,   serial  292162)  across  a  ten   ohm 

^^^il)  resistor  placed  in  series  with  the  circuit,  and 
dividing  the  voltage  by  ten  ohms.  The  resulting  graph  is 
shown  in  Figure  IV-7.  The  slope  of  the  line  gives  us  .- 

dl/dV  =  2.10  ±  0.005  A/V.  (IV-10) 

Combining  the  values  in  Equations  IV-9  and  IV-10  with  the 
previously  calculated  high  end  value  for  the  suspension 
stiffness  k  =  9.2±0.6X  10^  kg/s==  g  j  vcs  : 

Bl  =  k(dx/dV) (dV/dl)  =  21.5  ±  1.4  N/A.      (IV-11) 

This   15%  difference  between  the  values  of  Bl   using 

the  two  independent  methods  is  what  prompted  us  to  test   for 

Bl   and   k   using  the  static   loading   method.   Our   primary 

concern  is  that  we  would  expect  the  direct  measurements  of  B 
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and  1  to  give  us  the  highest  possible  value  for  Bl .  This  is 
because  the  maximum  magnetic  field  may  not  be  acting  on  the 
full  length  of  the  coil.  We  did  not  expect  the  value  of  Bl 
calculated  using  the  spring  constant  to  be  15%  higher  than 
the  direct  measurement  value. 

We  were  able  to  obtain  static  loading  values  for 
both  Bl  and  k  using  the  same  equipment  setup  and  the 
formulas : 

Bl  =  Amg/AI  and  k  =  -  Amg/Ax  (IV-12) 
For  these  measurements  the  voice  coil  is  driven  with  a  DC 
current.  We  then  add  previously  measured  masses  (putty  and 
soldering  wire)  to  the  voice  coil,  and  measure  the  resulting 
displacement  of  the  diaphragm  using  the  photonic  sensor.  To 
measure  Bl  we  determine  how  much  current  is  needed  to  bring 
the  voice  coil  back  to  its  null  displacement  (original 
output  voltage  reading  from  the  photonic  sensor)  after  the 
mass  is  added.  We  measure  k  by  keeping  the  current  constant 
and  measuring  the  difference  in  photonic  output  voltage 
(multiplied  by  the  calibration  constant  to  get  displacement) 
resulting  from  the  added  mass.  Table  IV-1  lists  the  results 
of  these  measurements . 

The  bracketed  value  in  each  k  and  Bl  pair  was 
calculated  using  the  measured  value  in  the  pair  and  the 
relationship : 

BlI  =  -kx.  (IV-13) 
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From  these  measurements  we  get  the  value  k  =  6.4  +  o.2  kg/s= 
as  reported  earlier,  and  the  transduction  coefficient: 

Bl  =  15.0  ±  0.2  N/A.  (IV-14) 

For  the  model  runs  we  used  a  value  of  Bl  =  19.0  for  the  high 
k  and  Bl  =  15.0  for  the  low  k. 

TABLE  IV-1.   STATIC  LOADING  OF  NIB  DRIVER  TO  DETERMINE 
SUSPENSION  STIFFNESS  AND  TRANSDUCTION  COEFFICIENT 

Am  (g)     AI  (mA) 
17.1 

17.1 

34.7       23.3 

34.7       22.3 

34.7 

(front  slope) 

3 .   Mechanical  Resistance 

The   next  parameter  we  measured  was   the   mechanical 
resistance  (R„)  using  the  relationship.- 

Rm  =  27rf«m«/Q.  (IV-15) 

We  attempted  to  find  the  driver  quality  factor  (Q)  using  the 

same   modified   impedance   technique   described   in   Section 

III.G.   As   can   be  seen  in  Figure  IV-8,   the  graph   of   the 

modified  impedance  magnitude  vs.   frequency  is  not  a   smooth 

resonance  curve,   as  it  was  for  the  290-16K  (see  Figure  III- 

9).   In   fact,   it   looks   like   the   superposition   of   two 
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different  impedance  curves,  the  first  one  shorter  and  fatter 
and  the  second  one  taller  and  thinner.  It  was  impossible  to 
accurately  measure  half -power  points  from  this  curve,  so  we 
chose  the  phase  angle  method.  (It  is  interesting  to  note 
that  the  JBL  2445J  (ceramic)  driver  did  have  a  smooth  modi- 
fied impedance  curve  similar  to  that  of  the  290-16K.)  Figure 
IV-9  shows  the  graph  of  phase  angle  vs .  frequency  over  a 
wide  range  of  frequencies.  Although  this  curve  generally 
approximates  an  arctangent,  the  segment  of  the  line  immedi- 
ately adjacent  to  the  resonant  frequency  (assumed  to  be  the 
frequency  where  0=0,  in  this  case  534.5  ±  0.5  Hz)  is  linear 
(for  a  frequency  range  of  about  ten  Hz) .  This  segment  is 
plotted  in  Figure  IV-10.  The  slope  of  the  line  is  dt>/df  =  - 
•42.0  ±  0.5  X  10~^s  (after  converting  degrees  to  radians).  As 
in  Section  II. G,  Q  is  given  by: 

0  =  -if^/2) (dO/df)  =  11.2  ±  0.2.  (IV-16) 
This  gives  a  value  for  R„  (using  the  moving  mass  in  air 
m^  =  7.3  ±  0.1  g): 

Rrn  =  2.34  ±  0.05  kg/s.  (IV-17) 

We  decided  that  a  more  accurate  measurement  of  Q  and 
Rm  could  be  found  by  making  the  Q  measurement  in  vacuum.  For 
this  we  used  the  open  circuit  free  decay  method.  This  method 
consists  of  sending  a  signal  to  the  driver  at  resonance  with 
an  HP  3314A  function  generator  (serial  A342044) ,  then 
simultaneously  switching  off  the  signal  and  triggering  a 
Nicolet  model  310  storage  oscilloscope  (serial  88004641)   to 
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record  the  free  decay.  A  sample  of  the  output  from  this 
process  is  shown  in  Figure  IV-11.  The  equation  for  a  simple 
harmonic  oscillator  experiencing  free  decay  is: 

A{t)  =  A„  exp  (-t/t)  sin(a>«t  +  -/)  ,  (IV-18) 
where  Ait)  is  the  amplitude  as  a  function  of  time,  A^  is  the 
initial  (maximum)  amplitude,  t  is  the  time  ( =  nT ,  or  the 
number  of  periods  times  the  period),  x  is  the  free  decay 
constant  (=  QT/jr,  see  Equation  III-29)  ,  ui„  is  the  angular 
resonant  frequency,  and  f  is  the  phase  angle.  We  recorded 
the  peak  amplitude  as  a  function  of  period  number,  and 
plotted  the  natural  logarithm  of  the  peak  amplitude  vs. 
period  number.  The  equation  representing  this  line  is: 

In  A{T)^..^  =  In  A«  -  nT/c.  {IV-19) 

The  slope  of  this  line  is  t/x  =  nT/^  =  nTr/Q.  Therefore  Q  can 
be  found  by  dividing  7t  by  the  slope  of  the  line.  Figure  IV- 
12  shows  a  semi-log  plot  for  the  NIB  driver  in  air.  The 
value  of  Q  for  this  plot  is  21.5,  which  is  almost  twice  the 
value  we  found  using  the  modified  impedance  phase  angle 
measurements  (Q  =  11.2).  Obviously  the  air  strongly  affects 
the  measurements .  The  scattering  of  the  data  points  in  this 
figure  is  probably  due  to  the  same  phenomena  that  causes  the 
double  peaks  in  the  impedance  curve.  Figure  IV-13  shows  the 
semi-log  plot  of  the  amplitude  vs.  period  number  for  vacuum. 
This  curve  is  much  smoother  and.  we  assume,  more  reliable. 
The  Q  for  this  data  is : 

0  =  7t/slope  =  19.8  ±  0.05.  (IV-20) 
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Figure  IV-ll 


Sample  of  Free  Decay  Output  from 
Nicolet  310  Storage  Oscilloscope 
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We  can  also  find  the  period  from  the  free  decay  curve  (and 
thus  the  resonant  frequency  as  the  inverse  of  the  period) . 
We  use  the  oscilloscope  to  measure  the  time  for  ten  cycles 
and  divide  this  number  by  ten.  For  the  vacuum  run  the 
resonant  frequency  was  839  ±  0.5  Hz.  These  values  (and  the 
moving  mass  =  3.2  g  in  vacuum)  give  us  the  mechanical  resis- 
tance in  vacuum: 

R,„  =  27tf^m.=  /Q  =  0.83  ±  0.01  kg/s .         (IV-21) 
This  is  the  value  we  used  for  the  model . 

Table    IV-2    gives   a   comparison   of    applicable 
parameters  for  the  three  drivers.   Notice  that  the  parameter 

Bl/R„«c:himo  is  almost  three  times  greater  for  the  NIB  driver 
than  the  290- 16K.  This  is  the  parameter  that  we  said  (Sub- 
section IV. A. 2)  we  wanted  to  maximize  to  get  the  highest 
thermoacoustic  efficiency.  Thus  the  driver  that  gives  us  the 
best  height  and  weight  characteristics  (see  photograph  of 
the  three  drivers  in  Figure  IV-14)  also  gives  us  the  best 
electroacoustic  efficiency. 

C.   MODELED  EFFICIENCY 

In  order  to  model  the  electroacoustic  efficiency  of  the 
driver,  we  need  to  know  what  acoustic  load  (impedance)  it 
will  see.  The  final  resonator  and  stack  composition  will  not 
have  been  determined  by  the  time  this  thesis  is  completed. 
We   therefore   assumed   a  range  of  possible  values   for   the 


110 


TABLE  IV-2.   COMPARISON  OF  PARAMETERS  FOR 
THE  THREE  DRIVERS 

Property  (units)       Altec  290      JBL  2445J       JBL  NIB 

Moving  mass,  m^  (kg)    0.0151         0.015  0.015 

Spring  ^constant  (N/m)   3. 52X10-*      -■'8.5X10^       -9.0X10^ 

Magnetic  Induction, 

B  (Tesla)  2.1  1.97  1.87 

Voice  Coil 

Diameter  (cm)         7.22  9.98  9.98 

Bl  (N/A)  21.5  18.5  17.5 

Voice  coil  DC 

resistance  (fi)        io.6  7.1  7.6 

Mechanical 

resistance  (kg/s)     2.77  --1.0  0.83 

Voice  coil 

inductance  (H)       1.0X10"^        1.14X10"*     1.14X10"^ 
(at  20  kHz) 

Mass  (kg)  14.5  13.8  3.1 

Height  (in)  6.0  5.5  2.5 

Magnetic  energy 

density  (J/cm^)       -^0.5  -^-0 . 5  3.2 

Bl/R„,«c=H,m«  510  1230  1410 

(m/Coul-kg) 
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Figure  IV-14 


Photographic  Comparison  of  the 
Three  Drivers — Altec  290-16K, 
JBL  2445J,  and  JBL  NIB. 
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acoustic  load  based  on  the  experimental  data  acquired  to 
date. 

The  model  parameters  that  relate  to  the  acoustic 
load  are  the  density  of  the  gas  (RHO) ,  the  speed  of  sound  in 
the  gas  (C) ,  the  resonator  quality  factor  (Q) ,  the  resonator 
resonant  frequency  (FO),  and  the  resonator  length  (L).  The 
model  was  designed  for  a  closed-ended  tude  of  uniform  cross- 
section.  The  STAR  resonator  has  an  open  tube  (with  a  bulb 
attached)  and  nonuniform  cross-section.  The  problem  of 
applying  the  existing  computer  model  to  the  STAR  resonator 
was  solved  by  scaling  the  STAR  parameters  to  make  it  look 
like  a  closed-ended  uniform  cross-section  tube. 
1 .   Input  Values  for  Pure  Helium 

For  pure  helium,  the  experimental  data  available  was 
the  range  of  resonator  quality  factor  values  (10-20),  and 
the  speed  of  sound  and  gas  density  at  room  temperature 
(about    300K)    and   ten   atmospheres   (Cr-t   =   1024    m/s , 

P''"'"  =  1.618  kg/m^)  .  Also  available  were  four  resonant  fre- 
quencies, two  for  a  short  resonator  at  warm  (580  Hz)  and 
cold  (520  Hz)  temperatures,  and  two  for  a  long  resonator  at 
warm  (520  Hz)  and  cold  (470  Hz)  temperatures.  Using  this 
information,  we  first  calculated  the  effective  lengths  of 
the  resonators  at  room  temperature  as  if  they  were  closed- 
ended,  uniform  cross-section  tubes.  For  this  we  used  the 
relationship : 

A  =  2L  =  c/f,   or   L.^^  =  c/2f,  (IV-22) 
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where  a  is  the  wavelength,  L  is  the  length  of  the  tube,  c  is 
the  speed  of  sound,  and  f  is  the  resonant  frequency.  This 
calculation  gave  us  effective  lengths  for  the  resonators: 
L,  =  0.883m  and  Li  =  0.985ni.  For  the  computer  runs  we  use  a 
convenient  value  of  L«^i-  =  0.90m.  We  then  calculated  the 
speed  of  sound  for  the  operating  temperature  of  approxi- 
mately 250K  using  the  relationship: 

c(T)  =  Cr-t  [T/Tr-t]  ^''=.  {IV-23) 

This  gives  us  a  speed  of  sound  c^-so  =  935  m/s .  If  we  use 
this  sound  speed  and  the  effective  resonator  lengths  to 
estimate  the  resonant  frequencies  for  250K,  we  get  f,  =  529 
Hz  and  f^  -  i\7A  Hz,  which  agree  with  the  measured  values.  We 
calculated  the  gas  density  at  the  operating  temperature 
using  the  relationship: 

p  =  ;^P/c=  =  1.94  kg/m=3,  (IV-24) 

where  z  is  the  ratio  of  the  specific  heats  (1.667  for 
helium)  and  P  is  the  pressure. 

Table  IV-3  lists  the  input  parameters  for  the  pure 
helium  computer  runs .  Note  that  for  the  computer  model .  the 
moving  mass  must  take  the  mass  of  the  reducer  cone  and 
bellows  into  account  (see  Figure  IV-15).  The  total  moving 
mass  is  expected  to  be  approximately  15  g. 
2.   Input  Values  for  Helium-Xenon 

There  was  somewhat  less  data  available  for  the 
helium-xenon  mixture  (12.5%  Xe) .  We  had  the  sound  speed  at 
300K   and   ten   atmospheres  (  Crt-  =  472   m/s).   We   also   had 
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TABLE  IV-3.  COMPUTER  MODEL  INPUT — 
NIB  DRIVER  WITH  HELIUM  AT  TEN  ATMOSPHERES 
AND  TEMPERATURE  =  250K 


Symbol  Meaning  Units       NIB 

RM  Mechanical  resistance 

M  Moving  mass 

SSUP  Spring  constant 

RHO  Density  of  gas 

C  Sound  speed 

S  Resonator  X-sec .  area 

BL  Driver  transduction  factor   N/A 

Q  Resonator  quality  factor 

FO  Resonator  resonance  f req . 

A  Back  volume  X-sec.  area 

VOL  Back  volume 

LE  Voice  coil  inductance 

RE  Voice  coil  DC  resistance 

I  Driver  current 

L  Resonator  length 


kg/s 

0.83 

kg 

0.015 

N/m 

6.6,  9.0  XIO-* 

kg/m^ 

1  .94 

m/s 

935.0 

m=2 

1.  14X10-^ 

N/A 

15.0,  19.0 

10. 15,20 

Hz 

519.4 

mz 

4.42X10-=' 

m=3 

3.28X10-^ 

H 

1  .  14X10--^ 

n 

7.6 

A 

1  .  0 

m 

0.90 
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(b) 


Figure  IV-15. 


Photograph  of  JBL  NIB  Driver 
with  (a)  Reducer  and  (b)  Bellows 
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resonant  frequencies  for  the  same  four  conditions  as  for  the 
pure  helium:  short-warm  (265  Hz),  short-cold  (230  Hz),  long- 
warm  (235  Hz),  and  long-cold  (215  Hz).  We  used  Equation 
IV-22  to  find  the  effective  resonator  lengths  (L,  =  0.391m 
and  Li  =  i.004m).  However,  once  again  we  used  a  convenient 
effective  length  of  0.90  m  for  the  computer  runs.  We  used 
Equation  IV-23  to  find  the  sound  speed  for  He-Xe  at  the 
operating  temperature,  c^-^o  =  431  m/s .  The  estimated  reso- 
nant frequencies  for  the  cold  temperature  were  f,=242  Hz  and 
f 1=215  Hz,  once  again  comparing  favorably  with  the  measured 
values.  To  find  the  density  of  the  mixture  at  the  operating 
temperature  we  used  a  form  of  the  ideal  gas  law: 

p  =  PM/RT,  (IV-25) 

where  M  is  the  atomic  mass  of  the  mixture,  R  is  the 
universal  gas  constant,  and  P  and  T  are  the  pressure  and 
temperature.  We  calculated  the  atomic  mass  of  the  mixture 
using  the  equation: 

Mm--x-  =  Mm.Xh,  +  Mx«x>s-,  (IV-26) 

where  the  M's  are  atomic  masses  (4  for  helium  and  131.3   for 

xenon)  and  the  x's  are  percentages  of  the  gases  in  the 
mixtures.  This  equation  gives  us  an  atomic  mass  for  the 
12.57o  Xe  mixture  of  19.9  g/mol ,  and  a  density  of  9.70  kg/m^ 
at  250K  and  ten  atmospheres.  Table  IV-4  lists  the  computer 
input  values  for  He-Xe,  including  a  range  of  Q's  from  5  to 
30.  (We  are  not  sure  what  to  expect  for  the  actual  value  of 
Q,  since  a  new  stack  will  also  be  used.) 
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TABLE  iy-4.  COMPUTER  MODEL  INPUT — 
NIB  DRIVER  WITH  HELIUM-XENON  AT  TEN  ATMOSPHERES 

AND  TEMPERATURE  =  250K 


Symbol  Meaning  Units       NIB 

RM  Mechanical  resistance 

M  Moving  mass 

SSUP  Spring  constant 

RHO  Density  of  gas 

C  Sound  speed 

S  Resor  -.or  X-sec .  area 

BL  Driver  trans .  factor 

Q  Resonator  quality  factor 

FO  Resonator  resonance  f req . 

A  Back  volume  X-sec.  area 

VOL  Back  volume 

LE  Voice  coil  inductance 

RE  Voice  coil  DC  resistance 

I  Driver  current 

L  Resonator  length 


kg/s 

0.33 

kg 

0.015 

N/m 

6  . 

6,  9.0  X10-* 

kg/m 

3 

9.  70 

m/s 

A31  .  0 

m= 

1 . 14X10-^ 

N/A 

15.0,  19.0 

5  , 

.  10 

, 15 , 20 , 25 , 30 

Hz 

239  .  4 

m= 

4  .  42X10-=* 

m^ 

3 . 23X10-^ 

H 

1  .  14X10--* 

il 

7.6 

A 

1  .0 

m 

0.90 
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3 .   Computer  Model  Results 

Figures  IV-16  through  IV-19  are  sample  output  plots 
for  pure  helium  with  Bl  =  15  N/A,  k  =  6.4  X10-*  kg/s^,  and 
Q  =  20.  Notice  that  the  peak  values  for  the  driver  and 
resonator  are  almost  matched  (the  driver  resonance  is  at  the 
lower  frequency) .  The  thermoacoustic  efficiency  for  this  set 
of  values  is  about  70.4%.  Figures  IV-20  through  IV-23  show 
the  corresponding  plots  for  the  helium-xenon  mixture.  This 
time  the  peaks  are  not  matched  (the  driver  peak  is  off  the 
scale)  and  the  electroacoustic  efficiency  is  about  16.6%. 
The  primary  reason  for  the  difference  between  the  two  gases 
is  the  lower  sound  speed,  and  therefore  lower  resonant  fre- 
quency, in  the  He-Xe  mixture. 

Table  IV-5  lists  the  maximum  value  of 
electroacoustic  efficiency  for  each  of  the  test  runs.  The 
consistently  low  efficiencies  for  the  He-Xe  emphasize  the 
importance  of  matching  the  two  resonant  peaks.  We  recommend 
that  the  driver  frequency  be  adjusted  to  match  the  resonator 
operating  frequency.  This  can  best  be  done  by  reducing  the 
stiffness  of  the  titanium  surround  (see  Section  III.H). 
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Figure  IV-16 


Computer  Model  Output  for  NIB  with 
Helium  (Bl  =  15,  Q  =  20)--Driver 
Voltage  vs.  Frequency  for  One  Amp. 
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Figure  IV-17 


Computer  Model  Output  for  NIB  with 
Helium  (Bl  =  15,  Q  =  20)--Piston 
Velocity  vs.  Frequency  for  One  Amp 
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Figure  IV-18 


Computer  Model  Output  for  NIB  with 
Helium  (Bl  =  15,  Q  =  20 ) --Acoustic 
Power  vs.  Frequency  for  One  Amp. 


122 


Figure  IV-19 


Computer  Model  Output  for  NIB  with  Helium 
(Bl  =  15,  Q  =  20)--Electroacoustic 
Efficiency  vs.  Frequency. 
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Figure  IV-20. 


Computer  Model  Output  for  NIB  with  Helium- 
Xenon  (Bl  =  15,  Q  =  20)--Driver  Voltage 
vs.  Frequency  for  One  Amp. 
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Figure  IV-21 . 


Computer  Model  Output  for  NIB  with  Helium- 
Xenon  (Bl  =  15,  Q  =  20)--Piston  Velocity 
vs.  Frequency  for  One  Amp. 
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Figure  IV-22 


Computer  Model  Output  for  NIB  with  Helium- 
Xenon  (Bl  =  15,  Q  =  20)--Acoustic  Power 
vs.  Frequency  for  One  Amp. 
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Figure  IV-23. 


Computer  Model  Output  for  NIB  with  Helium- 
Xenon  (Bl  =  15,  Q  =  20)--Electroacoustic 
Efficiency  vs.  Frequency. 
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TABLE  iy-5.   COMPUTER  MODEL  RESULTS  FOR  NIB  DRIVER-- 
PEAK  ELECTROACOUSTIC  EFFICIENCY  ( % ) 


Q. 

He 

(Bl=15) 

He 

(Bl=19) 

5 

10 

71  .3 

70.4 

15 

71.6 

74.3 

20 

70.4 

76.2 

25 

30 

5.1  5.9 

9.2  10.7 
13. 1  15.  1 
16.6  19.1 
19.6  22.7 
22.8  26.0 
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V.   STAR  DRIVER  CONSTRUCTION 

A.   INTRODUCTION 

This  chapter  describes  the  construction  of  the  STAR 
driver.  A  large  part  of  this  thesis  involved  the  design  of 
the  parts  necessary  to  integrate  the  NIB  driver  with  the 
resonator  and  GAS  cannister,  as  illustrated  in  Figure  1-4. 

The  major  design  work  was  for  the  driver  housing.  The 
requirements  for  the  driver  housing  are: 

1  -  lightweight, 

2  -  strong  (pressurized  to  210  psi), 

3  -  interface  to  GAS  canister  lid  (12  10-32  bolt  holes 

centered  at  9.5  inches  diameter), 

4  -  interface  with  resonator  (preferably  interchangeable 

with  earlier  test  resonators), 

5  -  helium  leak-tight  feed-through  ports  for  electrical 

leads  (driver  electric  power,  accelerometer  and 
microphone  outputs),  and 

6  -  hold  NIB  magnet  structure  securely. 

Figure  V-1  shows  a  cross-sectional  drawing  of  the  driver 
housing  (not  including  any  of  the  feed-through  ports,  which 
will  be  discussed  later).  Distinction  between  parts  is  shown 
by  different  hatching.  Part  A  is  the  NIB  magnet  structure 
which  was  discussed  in  Chapter  IV.  Part  B  is  the  main  body 
of   the   driver  housing,   machined  out  of  6061-T6   aluminum. 
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Figure  V-1 


Cross-sectional  Schematic  of  the  STAR  Driver 
Housing  Showing  (A)  NIB  Driver,  (B)  Driver 
Housing  Main  Body,  (C)  Aluminum  Pusher  Plate 
(D)  Pressure  Lid,  (E)  Leak  Test  Lid, 
(F)  Resonator  Vacuum  Housing,  and  (G)  Bolt 
Holes  for  Interface  with  Resonator. 
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Part  C  is  a  6061-T6  aluminum  pusher  plate  used  to  hold  the 
magnet  structure  securely  in  the  body  using  six  10/32  bolts. 
Part  D  is  the  pressure  lid,  sealed  a  l/16th  inch  diameter 
lead  o-ring  and  held  in  place  with  six  10/32  bolts.  Part  E 
represents  part  of  the  leak  test  lid  (see  Figure  C-1) .  When 
completed  this  will  be  replaced  by  the  GAS  can  lid  provided 
by  NASA.  The  driver  housing  must  match  up  with  the  GAS  can 
lid's  9.5  inch  diameter  10-32  bolt  circle.  Part  F  is  the 
resonator  section  housing.  The  arrows  labeled  G  point  to  the 
bolt  holes  which  form  the  driver-resonator  interface. 

B.   DESIGN  CONSIDERATIONS 

The  driver-resonator  system  will  be  pressurised  to  ten 
atmospheres  with  a  helium-xenon  mixture  (12.5%  xenon). 
NASA's  regulations  say  that  it  must  be  pressure  tested  to 
1.5  times  the  operating  pressure  (NASA,  1986a),  or  15  atmo- 
spheres (=1.52X10*  Pa  =  220  psi).  An  added  complication  of 
working  with  helium  (or  helium  mixtures)  is  that  it  diffuses 
through  almost  all  materials  except  metals.  As  a  result,  we 
are  not  able  to  use  rubber  o-rings  to  seal  metal  parts 
together,  but  must  use  lead  o-rings  instead. 

1 .   Electrical  Feed-Throughs 

NASA's  regulations  prohibit  us  from  leaking  gases 
into  space,  primarily  to  avoid  contamination  of  other  exper- 
iments and  for  the  safety  of  the  Shuttle  and  astronauts . 
This   regulation  is  important  not  only  for  the  major   o-ring 
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seals ,  but  also  for  the  microphone  and  accelerometer  outputs 
to  the  frequency  control  electronics  (see  Figure  1-5),  and 
connecting  the  audio  amplifier  to  the  driver.  There  are 
three  individual  wires  and  three  co-axial  connectors  that 
must  penetrate  the  driver  housing  without  allowing  leakage. 
The  other  housing  penetrations  are  a  DC  pressure  gage  for 
the  driver  housing  and  a  fill  line  from  a  helium-xenon 
bottle  which  will  stabilize  the  pressure  in  the  system 
against  temperature  changes  in  the  environment  and  those 
generated  by  the  refrigerator . 

Based  on  prior  experience  in  other  applications,  we 
decided  that  Emerson  &  Cuming's  STYCAST<tm>  2850FT  epoxy 
(Canton,  MA  02021)  would  be  the  best  choice  for  providing 
the  electrically  insulating  seal  for  the  wires  in  the 
housing.  STYCAST<-rM>  2850FT  is  a  quartz-filled  epoxy  formu- 
lation with  remarkably  good  over-all  general  properties.  In 
addition  to  having  excellent  electrical  grade  insulation 
properties,  good  adhesion,  and  excellent  resistance  to  chem- 
icals and  solvents,  2850FT  has  unusually  high  thermal  con- 
ductivity and  low  thermal  expansion.  This  combination  of 
properties  makes  this  epoxy  valuable  for  applications  where 
electrical  insulation  and  mechanical  protection  must  be 
maintained  while  coping  with  heat  transfer  considerations. 

One  definition  for  a  sealed  container  is  a  helium 
leak  rate  of  less  than  10--»  standard  cc's  per  second  (NASA. 
1986b) .   To   test   the   epoxy  we  designed  and  built   a   test 
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fixture  to  be  used  with  the  Veeco  MS-9  helium  leak  detector. 
A  cross-sectional  diagram  of  the  test  fixture  is  shown  in 
Figure  V-2 .  Part  A  is  a  copper  tube  with  1-1/8"  outer  diam- 
eter to  interface  with  the  leak  detector.  Soldered  to  the 
copper  tube  is  a  piece  of  solid  brass  stock  about  one  inch 
thick  (B)  with  a  0.052  inch  diameter  hole  bored  through  the 
center.  Epoxied  into  the  hole  is  a  0.025  inch  diameter  bare 
copper  wire  (C).  This  design  is  motivated  by  the  desire  to 
provide  as  long  a  travel  path  and  as  small  a  cross-section 
as  possible  for  the  helium  gas  molecules.  Part  D  is  a  hollow 
aluminum  housing  designed  to  hold  pressurized  helium.  We 
decided  that  a  rubber  o-ring  should  hold  the  helium  suffi- 
ciently well  for  a  short  term  leak  test.  Part  E  is  a  1/8" 
NPT  (National  Pipe  Thread)  hole  for  connecting  to  a  helium 
fill  line.  A  photograph  of  the  test  fixture  with  and  without 
the  housing  is  shown  in  Figure  V-3. 

The  copper  tube  and  brass  stock  were  soldered  toge- 
ther. We  boiled  the  assembly  to  remove  excess  particulates 
and  rinsed  it  with  alcohol.  We  then  etched  the  hole  in  the 
brass  and  the  copper  wire  with  nitric  acid  diluted  with 
distilled  water  (about  20%  acid)  to  provide  a  clean,  rough 
surface  for  the  epoxy .  Before  etching,  cellophane  tape  was 
placed  on  the  portion  of  the  brass  stock  where  the  o-ring 
would  be  sealing  (after  making  sure  that  the  tape  would  not 
dissolve  in  the  acid).  We  then  rinsed  the  assembly  with 
distilled  water  and  dried  it,   mixed  the  epoxy,  and  inserted 
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Figure  V-2 


Cross-sectional  Schematic  of  the  Helium 
Leak  Detector  Test  Fixture  Showing 

(A)  Copper  Tube  for  Leak  Detector  Interface, 

(B)  Solid  Brass  Stock,  (C)  Bare  Copper  Wire, 
(D)  Aluminum  Housing,  and  (E)  Connection  for 
Helium  Fill  Line. 
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Figure   V-3. 


Photograph  of  the  Helium  Leak  Detector 
Test  Fixture. 
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the  wire  into  the  hole  in  the  brass  stock,   making  sure   the 
epoxy  filled  the  gap. 

The  following  procedure  was  used  to  run  the  leak 
detector  tests:  first  zero  the  meter,  then  do  a  test  run  on 
the  sensitivity  calibrator.  The  sensitivity  calibrator  used 
was  a  Veeco  type  SC-4  serial  number  3179  with  a  leak  rate  of 

6.5X10~®  std  cc/sec  (standard  cubic  centimeters  per  second). 
We  would  then  remove  the  calibrator,  recheck  the  meter  for 
zero,  attach  the  test  fixture,  read  the  meter,  remove  the 
fixture,  and  recheck  for  zero. 

The  test  fixture  was  initially  pressurized  to  115 
psi  with  helium.  We  tried  to  test  it  immediately  but  found 
that  the  detector  would  not  draw  enough  of  a  vacuum  to  check 
for  helium.  We  thought  that  there  might  be  outgassing  from 
flu.x  or  something  in  the  tube.  So  we  pulled  a  vacuum  on  the 
tube  overnight  and  the  next  day  it  worked  fine.  We  tested 
the  fixture  on  the  first,  third,  and  sixth  days  after  it  had 
been  pressurized.  The  numerical  results  are  listed  in  Table 
V-1.  Even  after  we  had  waited  six  days  for  the  helium  to 
diffuse  through  the  epoxy,  the  fixture  leakage  rate  was  an 
order  of  magnitude  less  than  the  calibrator.  The  leakage 
rate  was  therefore  four  orders  of  magnitude  less  than  that 
allowed  by  NASA  to  be  considered  leak-tight  (1X10--*  std 
cc/sec).  After  the  third  leak  test  we  checked  the  pressure 
in  the  aluminum  chamber  and  found  it  to  be  within  a  few 
percent   of   the  original  pressure.   We  therefore   concluded 


136 


1 . 1 
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0.08 

1.  1 

0.  10 

that  the  bare  wire-epoxy  method  with  a  length/diameter  ratio 
of  about  20  could  be  used  for  the  driver  housing  penetrating 

electrical  feed-throughs . 

TABLE  V-1.  HELIUM  LEAK  DETECTOR  TEST  RESULTS 
ON  COPPER  WIRE-EPOXY  FIXTURE. 

Number  days  after  null  meter   calibrator   copper  wire  feed 

pressurization  reading      reading     -through  reading 

1  0 

3  0 

6  0 

Figure  V-4  shows  a  close-up  of  the  feed-through  plug 
for  the  wire  leads.  It  will  house  a  total  of  five  wires 
{with  space  available  for  a  sixth  if  necessary).  Two  of  the 
leads  will  be  for  the  voltage  supply  to  the  voice  coil . 
These  leads  will  consist  of  bare  copper  wires  epoxied  into 
the  small  holes  seen  in  the  figure.  Each  wire  will  then  be 
soldered  to  a  6/32  phillip-head  screw  with  a  1-72  blind- 
tapped  hole  in  it.  The  screw  will  be  glued  into  an  8/32 
hole,  forming  a  mechanically  stable,  demountable,  electrical 
terminal.  The  accelerometer  and  microphone  leads  will  also 
have  bare  copper  wire  epoxied  into  the  small  hole.  and  will 
then   be   connected  to  Microdot 'tm>  coaxial  fittings   (Malco 
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Figure   V-^^ 


Driver    Housing    Feed-Through    Plug 
for    the    Electrical    Leads. 
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part   number   031-0001-0001,   South   Passadena ,   CA   91030). 
Figure  V-5  shows  an  example  of  each  of  these  leads . 
2 .   Measurement  Devices 

The  microphone  and  accelerometer  used  by  STAR  are 
similar  to  those  used  by  Hofler  (1986  and  1988).  In  fact, 
the  resonator  end  of  the  driver  housing  was  designed  to 
accept  Hofler 's  resonator,  calibration  volume,  and  instru- 
mentation for  testing.  The  microphone  used  to  measure  the 
acoustic  pressure  in  the  resonator  is  a  small  quartz-crystal 
( Valpey-Fisher ,  Hopkinton,  MA  01748)  or  PZT  dynamic  pressure 
transducer  (Piezoelectric  Products  lot  number  4-9-8A, 
F013184,  Metuchen,  N J ) .  It  is  mounted  in  the  driver  housing 
where  the  housing  connects  to  the  resonator  flange.  The 
microphone  will  be  calibrated  by  comparison  to  a  standard 
transducer,  which  will  have  been  previously  calibrated  via  a 
reciprocity  technique.  The  microphone  will  be  connected  to  a 
miniature  FET  follower  circuit.  The  FET  buffer  (Eletec  model 
304  impedance  converter,  Daytona  Beach,  FL  32020)  converts 
the  microphone's  high  impedance  output  signal  to  a  low 
impedance  to  reduce  sensitivity  to  electrical  noise  and 
maintain  the  output  voltage  and  phase  while  driving  the 
cables  necessary  to  reach  the  signal  processing  electronics. 
The  buffer's  DC  power  supply  and  output  cables  are  complete- 
ly shielded  through  the  driver  housing  to  the  buffer  (see 
Figures  C-3  and  C-7  for  location  of  microphone  and  FET  in 
driver  housing). 
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(a) 


(b) 


Figure  V-5. 


Example  of  Wire  Lead  Connectors 
for  Feed-Through  Plug:  (a)  Single 
Wire  and  (b)  Co-axial  Cable. 
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The  Endevco  model  22  Picomin ''^'^'  acceierometer  is 
attached   to   the  side  of  the  bellows  stiffener   facing   the 

driver.  The  integrated  acceierometer  signal  will  be  used  as 
a  measure  of  the  bellows  velocity,  which  will  be  calibrated 
by  driving  a  small  cavity  of  calculable  acoustic  impedance. 
Thus  the  product  of  the  outputs  from  the  microphone  and 
integrated  acceierometer  can  be  combined  to  yield  the  acous- 
tic power  and  their  relative  phase  can  be  used  to  keep  the 
driver  frequency  at  the  tube  resonance.  The  operation  of  the 
refrigerator  at  resonance  is  essential  both  to  optimize  the 
acoustic  power  delivered  to  the  "stack"  and  to  guarantee 
that  the  "stack"  is  always  located  at  the  proper  position  in 
the  standing  wave  as  discussed  in  Sub-section  II. A. 4. 

The  DC  pressure  transducer  is  an  Omega  ^'^'''^  PX-30 
designed  for  high  pressure  applications  requiring  long-term 
stability.  Figure  V-6  shows  a  drawing  of  the  pressure  trans- 
ducer and  its  driver  housing  fitting. 

The   helium   port   in  the  housing  is   provided   to 
empty  and  fill  the  driver  and  resonator  and  so  that  a  stable 
pressure   may   be   maintained  in  the   driver   and   resonator 
through   the  ballast  volume  provided  by  the   reservoir   (see 
Figure  1-4 ) . 

There  is  a  1/32  inch  diameter  hole  filled  with  a  ten 
thousandths  inch  inner  diameter  half  inch  long  capillary  in 
the  driver  housing  where  it  connects  to  the  resonator.  The 
function   of   the  capillary  is  to   equilibriate   the   static 
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pressure  between  the  driver  and  resonator  during  filling  and 
non-operating  times  while  maintaining  a  pressure  seal  at 
acoustic   operating  frequencies. 

A  complete  set  of  dimensioned  drawings  for  the 
driver  housing  and  its  auxiliary  equipment  is  given  in 
Appendix  C.  A  list  of  materials  for  the  driver  portion  of 
the  STAR  is  given  in  Appendix  D. 
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VI.   CONCLUSION  AND  RECOMMENDATIONS 

The  objective  of  the  STAR  project  is  to  test  and  space 
qualify  a  new  continuous  cycle  cryogenic  refrigeration 
system  for  cooling  of  sensors  and  electronics  which  is  based 
upon  the  newly-discovered  thermoacoustic  heat  pumping 
effect.  The  new  refrigerator  has  no  sliding  seals,  a  cycle 
frequency  of  about  300  Hz,  and  uses  acoustic  resonance  to 
enhance  overall  efficiency.  This  thesis  is  concerned  speci- 
fically with  the  design  and  testing  of  the  electrodynamic 
transducer  which  is  responsible  for  the  electro-acoustic 
power  conversion.  A  computer  model  of  the  driver/resonatoi- 
system  was  presented  along  with  the  techniques  for  measure- 
ment of  the  electrical  and  mechanical  parameters  used  as 
input  for  the  model . 

A  final  driver  design  (including  dimensional  drawings) 
utilizing  a  modified  JBL  2'450J  neodymium-iron-boron 
compression  driver  and  associated  leak-tight  electrical 
f eed-throughs ,  microphone,  accelerometer ,  pressure  gage, 
pressure  housing,  and  resonator  interface  is  provided. 

Since  we  only  received  the  NIB  driver  in  late  December 
1987,  there  wasn't  time  to  perform  enough  parameter 
measurements  to  obtain  accurate  values  for  the  suspension 
stiffness   (k)  and  transduction  coefficient  (Bl).   The   next 
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group  of  students  should  find  these  values  for  use  in  the 
model  and  efficiency  calculations. 

The  resonator  experiments  (not  reported  here)  indicate 
that  Helium-Xenon  may  be  the  preferred  gas  mixture  for  the 
STAR.  If  this  is  the  case,  the  NIB  driver  resonant  frequency 
will  have  to  be  lowered  for  it  to  match  the  operating 
(resonant)  frequency  of  the  resonator.  This  could  be 
accomplished  by  adding  moving  mass  to  the  voice  coil- 
reducer-bellows  assembly.  but  that  would  decrease  the 
velocity  and  lower  the  electroacoustic  efficiency.  The 
preferred  option  for  driver  resonance  adjustment  is  to 
decrease  the  suspension  stiffness  by  cutting  away  part  of 
the  surround  (see  Section  III.H).  Since  the  NIB  voice  coil 
has  a  titanium  surround  and  diaphragm  this  operation  will 
have  to  be  performed  carefully.  The  titanium  diaphragm 
(dome)  will  also  have  to  be  removed.  (The  dome's  function  is 
performed  by  the  reducer  and  bellows . ) 

The  task  for  the  next  group  of  students  is  to  assemble 
the  NIB  driver  with  its  housing  and  accessories  (feed- 
throughs ,  microphone,  accelerometer ,  etc.)  and  perform 
efficiency  measurements  on  the  ground  similar  to  those  that 
will  be  performed  in  space.  This  series  of  measurements  will 
closely  resemble  those  performed  on  the  Altec  290-16K  driver 
as  described  in  Chapter  III. 
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APPENDIX  A.   COMPUTER  MODEL  CALCULATION  PROGRAM 


PROGRAM   FRIDGE 

C  INTRODUCTION: 
C 

C  THE  PURPOSE  OF  THIS  PROGRAM  WRITTEN  BY  LT .  MICHAEL  P. 
C  SUSALLA  IS  TO  CALCULATE  THE  DATA  POINTS  AND  PLOT  THE 
C  PERFORMANCE  CHARACTERISTICS  OF  A  CONCEPTUAL  MODEL  OF  A 
C  THERMOACOUSTIC  HEAT  ENGINE.  THE  COMMENTS  INCORPORATED  INTO 
C  THIS  PROGRAM  AND  THE  COMMENTS  INCORPORATED  IN  THE  PLOTTING 
C  PROGRAM  "GRAFRI  FORTRAN  Al",  WILL  ALLOW  THE  USER  TO  PRO- 
C  DUCE  FOUR  GRAPHICAL  PLOTS  WHICH  WILL  CLOSELY  APPROXIMATE 
C  THE  PERFORMANCE  OF  AN  ACTUAL  THERMOACOUSTIC  HEAT  ENGINE. 
C  THE  GOAL  OF  THE  PROGRAM/SIMULATIO^^  ^S  TO  ALLOW  THE 
C  OPERATOR  TO  VARY  DIFFERENT  PARAMETERS  TO  ACHIEVE  WHAT  IS 
C  BELIEVED  TO  BE  AN  OPTIMUM  CONDITION  WHERE  THE  ENGINES 
C  DRIVER  RESONANCE  IS  MATCHED  WITH  THE  TUBE  RESONANCE. 

C  GENERAL  INSTRUCTIONS: 

C 

C      THIS  PROGRAM  IS  SET  UP  TO  PROVIDE  FOUR  PLOTS:   VOLTAGE 

C  VERSUS  FREQUENCY,  PISTON  VELOCITY  VERSUS  FREQUENCY.  ACOUS- 

C  TIC  POWER  VERSUS  FREQUENCY,  AND  ELECTROACOUSTIC  EFFICIENCY 

C  VERSUS   FREQUENCY.   IN  GENERAL  THE  PROGRAM  HAS  HARD   INPUT 

C  DATA   SUCH   THAT  ALL  YOU  MUST  DO  TO  GET  A  SET  OF  PLOTS   IS 

C  COMPILE  "GRAFRI  FORTRAN  Al"  IN  FORTVS  AND  THEN  RUN  "FRIDGE 

C  FORTRAN  Al"  IN  WF77  FROM  FLIST  ON  THE  NPS  MAINFRAME   USING 

C  A  TEKTRONICS  618  TERMINAL.  THIS  ACTION  WILL  GIVE  YOU  A  SET 

C  OF  PLOTS  FOR  THE  INPUT  DATA.  TO  TRY  OTHER  DATA  YOU  MUST  GO 

C  TO   XEDIT  ON  THE  "FRIDGE  FORTRAN  Al"  PROGRAM  AND  VARY   THE 

C  INPUT  DATA.   AS  LONG  AS  YOUR  CHANGES  ARE  RESTRICTED  TO  THE 

C  "FRIDGE"   PROGRAM  THE  ONLY  ACTION  NECESSARY  TO  PRODUCE  NEW 

C  PLOTS  IS  TO  RUN  THE  PROGRAM  IN  WF77 .  IF  YOU  FIND  IT  NECES- 

C  SARY  TO  CHANGE  THE  PLOTTING  PROGRAM,   "GRAFRI".   FOR  WHAT- 

C  EVER  REASON,   IT  MUST  BE  RECOMPILED  IN  FORTVS  PRIOR  TO  THE 

C  NEXT   ATTEMPT  TO  PLOT  OR  THE  CHANGES  WILL  NOT  BE   SEEN   BY 

C  THE  MACHINE. 

C 

C  INPUT  SECTION: 

C      A.   DEFINITION  OF  VARIABLES: 

C  ALL  OF  THE  VARIABLES  HAVE  BEEN  ESTABLISHED  AS  REAL 

C  OR   COMPLEX   QUANTITIES  BASED  ON  THE  REQUIREMENTS   OF   THE 

C  ARITHMETIC   ALGORITHOM  OF  THE  CONCEPTUAL  MODEL.     EACH   OF 

C  THE  VARIABLES  WILL  NOW  BE  DEFINED. 

C      LCOIL=THE  IMAGINARY  PART  OF  THE  IMPEDANCE  OF  THE  VOICE 

C  COIL  OF  THE  DRIVER,  AND  IS  CALCULATED  TO  BE  W*LE. 

C      K=THE   PROPAGATION   CONSTANT  OF  THE  WAVES  IN  THE   TUBE 

C  AND  IS  GIVEN  BY  W/C . 
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C      IMZM=THE  IMAGINARY  PART  OF  THE  MECHANICAL  IMPEDANCE  OF 

C  THE  TUBE/DRIVER  COMBINATION.   IT  IS  GIVEN  BY  A  RATHER  LONG 

C  AND   COMPLEX   EQUATION  IN  THE  CALCULATION  SECTION   OF   THE 

C  PROGRAM. 

C      MAGU=THE   MAGNITUDE  OF  THE  COMPLEX   QUANTITY   "U" ,   IN 

C  METERS/SEC. 

C      IMZT=THE   IMAGINARY  PART  OF  THE  TUBE  MECHANICAL   IMPE- 

C  DANCE  AND  IS  GIVEN  BY  SUBTRACTING  THE  DRIVER  PORTIONS  FROM 

C  IMZM. 

C      M=THE  MOVING  MASS  OF  THE  SYSTEM,  AND  IS  THE  SUM  OF  THE 

C  VOICE   COIL   AND  THE  DIAPHRAGM  MASSES  OF   THE   DRIVER,   IN 

C  KILOGRAMS. 

C      L=THE  LENGTH  OF  THE  TUBE  IN  METERS. 

C      LE=THE  INDUCTANCE  OF  THE  VOICE  COIL,  IN  HENRYS. 

C      I=THE  INPUT  CURRENT  FROM  THE  SOURCE,  IN  AMPS. 

C      MAGV  =  THE   MAGNITUDE   OF  THE  COMPLEX  QUANTITY   "V,   IN 

C  VOLTS. 

C      ZM=THE  MECHANICAL  IMPEDANCE  OF  THE  TUBE/DRIVER  SYSTEM. 

C  AND  IT  IS  MADE  UP  OF  THE  REAL  PART,   REZM,   AND  THE  IMAGI- 

C  NARY  PART.  IMZM. 

C      ZE=THE   ELECTRICAL  EQUIVALENT  OF  THE  MECHANICAL   IMPE- 

C  DANCE  OF  THE  TUBE/DRIVER  SYSTEM.   IT  IS  GIVEN  BY  TRANSFOR- 

C  MING   ZM  ACROSS  A  TRANSFORMER  USING  THE  BL  COEFFICIENT   OF 

C  THE  DRIVER. 

C      ZCOIL=THE   COMPLEX  ELECTRICAL  IMPEDANCE  OF   THE   VOICE 

C  COIL. 

C      ZTOT=THE   TOTAL   ELECTRICAL  IMPEDANCE  OF   THE   SYSTEM, 

C  COMPLEX. 

C      V=THE   COMPLEX   VOLTAGE   OF  THE  SOURCE   NECESSSARY   TO 

C  PRODUCE  THE  SELECTED  INPUT  CURRENT. 

C      U=THE   VELOCITY   OF   THE   PISTON ,( DIAPHRAGM ) ,   OF   THE 

C  DRIVER  IN  M/S. 

C      ZTUBE=THE   COMPLEX   MECHANICAL  IMPEDANCE  OF   THE   TUBE 

C  ITSELF. 

C      P=THE  COMPLEX  PRESSURE  OF  THE  SYSTEM,  IN  NEWTONS/M* *2 . 

C      PI=A   CONSTANT   FOR  CONVERSIONS  TO   AND   FROM   ANGULAR 

C  UNITS. 

C      RM=THE  MECHANICAL  RESISTANCE  OF  THE  MOVING  MASS  OF  THE 

C  SYSTEM. 

C      SSUP=THE  MECHANICAL  STIFFNESS  OF  THE  SYSTEM  DUE  TO  THE 

C  DRIVER  SUSPENSION. 

C      RHO=THE   DENSITY  OF  THE  GAS  IN  THE  TUBE  AT  A  SPECIFIED 

C  MEAN  TEMPERATURE. 

C      C=THE   SPEED   OF   SOUND   FOR   THE   SPECIFIED   GAS   AND 

C  TEMPERATURE. 

C      S=THE  CROSS  SECTIONAL  AREA  OF  THE  TUBE. 

C      BL=A   COEFFICIENT  OF  THE  DRIVER  THAT  IS  BASED   ON   THE 

C  STEADY   AIR  GAP  FLUX  DENSITY  OF  THE  DRIVER  AND  THE   LENGTH 

C  OF  ITS  VOICE  COIL. 

C      Q=THE  QUALITY  FACTOR  OF  THE  TUBE. 

C      FO=THE  RESONANCE  FREQUENCY  OF  THE  TUBE. 

C      A=THE   CROSS  SECTIONAL  AREA  OF  THE  BACK  VOLUME  OF   THE 
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C  DRIVER. 

C      VOL=THE  VOLUME  OF  THE  BACK  VOLUME  OF  THE  DRIVER. 

C      RG=THE   INTERNAL   ELECTRICAL  RESISTANCE  OF   THE   POWER 

C  SOURCE. 

C      RE=THE  ELECTRICAL  RESISTANCE  OF  THE  VOICE  COIL. 

C      SVOL=THE   MECHANICAL   STIFFNESS  OF  THE  SYSTEM  DUE   THE 

C  BACK  VOLUME  OF  THE  DRIVER. 

C      WO=THE   RESONANCE   FREQUENCY  OF  THE   TUBE   IN   ANGULAR 

C  UNITS,  RAD/SEC. 

C      ALPHA=A   CHARACTERISTIC  QUANTITY  OF  THE  TUBE  THAT  IS  A 

C  FUNCTION  OF  THE  RESONANT  FREQUENCY,  THE  Q  AND  THE  SPEED  OF 

C  SOUND. 

C      F=FREQUENCY.  IN  HERTZ. 

C      W=THE  FREQUENCY  CONVERTED  TO  RAD/SEC. 

C      FORCE=THE   FORCE  ON  THE  MOVING  MASS  DUE  TO  THE   INTER- 

C  ACTION   OF   THE  VOICE  COIL  MAGNETIC  FIELD  AND   THE   DRIVER 

C  MAGNET. 

C      REZT=THE   REAL   PART  07  THE   COMPLEX   MECHANICAL   TUBE 

C  IMPEDANCE. 

C      PU^THE   PRESSURE  TIMES  THE  VOLUME   VELOCITY,   ACOUSTIC 

C  POWER. 

C      PWR=THE  ELECTRICAL  POWER  NECESSARY  TO  DRIVE  THE  SYSTEM 

C  AT  THE  SPECIFIED  INPUT  CURRENT. 

C      PUPWR=A  POWER  RATIO,  ACOUSTIC  TO  ELECTRICAL.   AN  SFFI- 

C  CIENCY. 

C 

C      B.DATA  DESCRIPTION. 

C  FOR  THE  SELECTED  INPUT  DATA  THAT  IS  HARD  INPUT  THE 

C  FOLLOWING  VALUES  WERE  EMPIRICALLY  DETERMINED:  M,  L,  Q,  FO , 

C  A,   VOL,   BL,  S,  SSUP,  AND  RE.  RHO  AND  C  WERE  SELECTED  FOR 

C  HELIUM  GAS  AT  10  ATMOSPHERES  AND  A  MEAN  TEMPERATURE  OF  300 

C  DEGREES   KELVIN.    RM   WAS  CALCULATED  FROM   A   THEORETICAL 

C  EQUATION  USING  RHO , C  S  AND  THE  DRIVER  RESONANT   FREQUENCY. 

C  RG  WAS  SET  TO  ZERO  TO  MODEL  AN  IDEAL  POWER  SOURCE.   LE  WAS 

C  SET  TO  ZERO  FOR  EASE  OF  CALCULATION,  BUT  WAS  TESTED  DURING 

C  RUNNING   TO   VERIFY  THAT  IT  HAD  LITTLE  IF  ANY  EFFECT   WHEN 

C  INCLUDED.   A  CURRENT,   I,   OF  ONE  AMPERE  WAS  SELECTED   FOR 

C  EASE  OF  ANALYSIS. 

C      SINCE   YOU  MAY  WANT  TO  CHANGE  MY  INPUT  DATA  I  WILL  PUT 

C  A   COPY   IN  THE  COMMENTS  SO  THAT  YOU  CAN  RETURN  TO   IT   IF 

C  DESIRED.   MY  ANALYSIS  HAS  SHOWN  THE  BEST  BACK  VOLUME  TO  BE 

C  9.50E-8,  AND  THE  BEST  RM  TO  BE  1.535  KG/S.   THIS  VALUE  FOR 

C  RM  WAS  DETERMINED  AFTER  REVIEW  OF  PLOTS  SHOWED  HIGHER  THAN 

C  EXPECTED  VOLT.AGES  AND  IT  WAS  DETERMINED  BY  MEASUREMENT   IN 

C  THE  LABORATORY. 

C 

C 

C      DATA  RM,M,SSUP,RH0,C/3. 53E-3,0 .0045, 2.6E4. 1 .68, 1000/ 

C      DATA  S, BL,Q,FO, A/1 .521E-3, 8, 15,800,7. 13E-5/ 

C      DATA  VOL.RG,LE,RE, I ,L/2.715E-6,0 ,0,8, 1 ,0.625/ 
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REAL  LCOIL , K , IMZM , MAGU , IMZT , M , L , L£ , I , MAGV 

COMPLEX  ZM , ZE , ZCOI L , ZTOT , V , U . ZTUBE , P 

PI=3. 141592654 

DATA  RM,M,SSUP,RH0.C/3.53E-3, 0 .0045, 2. 6E4, 1 . 63 , 1000/ 

DATA  S,BL,Q,F0,A/1 .521E-3,8, 15,800,7. 13E-5/ 

DATA  VOL , RG , LE , RE , I , L/2 . 715E-6 .0,0,3,1,0. 625/ 

SVOL={A*A*RHO*C*C)/VOL 

W0=2*PI*F0 

ALPHA=W0/(2*Q*C) 


C  CALCULATION  SECTION: 

C  THE  NEXT  SIX  LINES  CREATE  AN  OUTPUT  FILE  FOR  THE  DATA 
C  THAT  IS  NECESSARY  TO  PRODUCE  THE  DESIRED  PLOTS,  AND  PLACES 
C  HEADINGS  ON  THE  DATA  AS  APPROPRIATE. 


0PEN{UNIT=69,FILE='FRIDGE  OUTPUT  Al ' , IOSTAT=IOS ) 
WRITE(69,*)'     MASS         VOLUME' 
WRITE(69,  DM,  VOL 

WRITE(69,*)'    FREQUENCY      PIST  VEL   PRESS*PIST  VEL   PWR  RATIO 
*     VOLTAGE' 
1  F0RMAT(E13.4,E13.4) 


C  THE  FOLLOWING  DO  LOOP  PROVIDES  FOR  THE  CALCULATION  OF 
C  THREE  HUNDRED  DATA  POINTS  FOR  EACH  OF  THE  DESIRED  PLOTS, 
C  USING  A  FREQUENCY  RANGE  OF  ONE  QUARTER  TUBE  RESONANCE  TO 
C  THREE  HALVES  TUBE  RESONANCE. 


DO  10  J=l,300 

F=(F0/4)+{ ( J-1)*(FO/240) ) 


C  THE  NEXT  THREE  LINES  CALCULATE  VARIOUS  QUANTITIES. 
C  EXPLAINED  ABOVE.  THAT  ARE  FUNCTIONS  OF  FREQUENCY,  FOR  USE 
C  IN  FUTURE  CALCULATIONS. 


W=2*PI*F 
K  =  W/C 
LCOIL=W*LE 


C        THE   NEXT  NINE  LINES  CALCULATE  THE  REAL  AND   IMAGINARY 
C  PARTS   OF  THE  TUBE/DRIVER  MECHANICAL  IMPEDANCE  AND  COMBINE 
C  THEM  INTO  A  COMPLEX  QUANTITY,  ZM . 


REZM  =  RM+(  (  (RHO*C*S* ALPHA *COS(K*L)  "*SIN(K*L)  )  +  (K*RHO*C*S 
*SINH( ALPHA* L) *COSH ( ALPHA  * L )  )  )/(K* (1+(ALPHA/K) **2)  * (  (  (  ( 
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*  SIN(K*L) ) **2) * { ( COSH ( ALPHA *L) ) **2) )+( ( (COS(K*L) ) **2) * ( ( 

*  SINH{ALPHA*L) ) **2) ) ) ) ) 

IMZM= { W  *  M ) - ( SSUP/W ) - ( SVOL /W )  +  (  (  { RHO  *  C  *  S  *  ALPHA  *  SINH { ALPHA  * 

*  *COSH(ALPHA*L)  )  -  (  K*RHO*C*S*COS  (  K*L  )  *SIN(K*L)  ))/(!<*  (1+(ALP 

*  /K) **2) * ( ( ( (SIN{K*L) ) **2) *( ( COSH ( ALPHA* L) ) ^*2) )+( ( (COS(K* 

*  **2) * ( (SINH(ALPHA*L) ) **2) ) ) ) ) 
ZM=CMPLX ( REZM , IMZM ) 


C  THE  NEXT  THREE  LINES  CONVERT  ZM  TO  AN  ELECTRICAL 
C  QUANTITY  AND  THEN  COMBINE  ALL  OF  THE  ELECTRICAL  IMPEDANCES 
C  TO  GET  A  TOTAL  IMPEDANCE. 


ZE=(BL**2)/ZM 
ZCOIL=CMPLX ( RE , LCOIL ) 
ZTOT=RG+ZCOIL+ZE 


C  THE  NEXT  SIXTEEN  LINES  CALCULATE  THE  VALUES  TO  BE 
C  PLOTTED,  AND  PUT  THEM  IN  THE  OUTPUT  FILE,  AND  THEN  CLOSE 
C  THE  FILE. 


V=I*ZTOT 
FORCE=BL*I 
U=FORCE/ZM 
MAGU=CABS(U) 
REZT=REZM-RM 

IMZT=IMZM- ( W*M ) + ( SSUP/W ) + ( SVOL/W ) 
ZTUBE=CMPLX ( REZT , IMZT ) 
P=(ZTUBE*U)/S 

PU=0.5* (REAL(P*CONJG(U) ) *S) 
PWR=0.5*(REAL(I*CONJG(V) ) ) 
PUPWR=PU/PWR 
-  MAGV=CABS(V) 

WRITE ( 69 , 100 ) F , MAGU , PU , PUPWR , MAGV 
100     FORMAT ( El 3. 5, El 3. 5, El 3. 5, E13. 5, El 3. 5) 
10  CONTINUE 
CLOSE (69) 


C  PLOTTING  SECTION: 

C      THE  NEXT  TWO  LINES  INTERFACE  WITH  MY  PLOTTING  PROGRAM, 

C  "GRAFRI  FROTRAN  Al",   AND  THE  DISSPLA  ROUTINE  TO  PLOT   THE 

C  DATA. 


PAUSE  'PRESS  ENTER  TO  START  PLOTTING  ROUTINE' 
CALL  EXCMS('EXEC  DISSPLA  GRAFRI') 
STOP 
END 
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APPENDIX  B.   COMPUTER  MODEL  GRAPH  PROGRAM 


PROGRAM  GRAFRI 

C  THIS  IS  LT.  MICHAEL  P.  SUSALLA * S  PLOTTING  PROGRAM  FOR 
C  USE  IN  CONJUNCTION  WITH  "FRIDGE  FORTRAN  Al"  IN  "THE  THER- 
C  MOACOUSTIC  HEAT  ENGINE  SIMULALTION" .  THIS  PROGRAM  MUST  BE 
C  COMPILED  IN  FORTVS  PRIOR  TO  RUNNING  "FRIDGE  FORTRAN  Al"  IN 
C  WF77.THIS  PROGRAM  IS  FOR  USE  IN  DISSPLA  AND  CALLS  VARIOUS 
C  DISSPLA  SUBROUTINES.  THESE  VARIOUS  SUBROUTINES  ARE  USED 
C  MORE  THAN  ONCE  IN  THE  PROGRAM  BUT  WILL  ONLY  BE  DESCRIBED 
C  ONCE  EACH.  IN  ORDER  TO  USE  DISSPLA  YOU  MUST  HAVE  ARRAYS 
C  OF  DATA  POINTS  TO  BE  PLOTTED.  MY  ARRAYS  ARE  ALL  DIMEN- 
C  SIGNED  TO  300  POINTS  SINCE  THAT  IS  THE  NUMBER  CALCULATED 
C  BY  "FRIDGE  FORTRAN  Al".  YOU  SHOULD  BE  AWARE  THAT  IF  YOU 
C  CHANGE  THE  NUMBER  OF  POINTS  TO  BE  CALCULATED,  AND  THE 
C  NUMBER  IS  GREATER  THAN  300,  YOU  MUST  REDIMENSION  THE 
C  ARRAYS.  TWO  CHARACTER  STRINGS  ARE  USED  IN  THIS  PROGRAM. 
C  NEW1$  IS  USED  IN  RENAMING  FILES  TO  THE  FORMAT  NECESSARY 
C  FOR  DISSPLA,  AND  DUMMY$  IS  USED  TO  READ  THE  HEADINGS  OUT 
C  OF  THE  DATA  FILE.  ALL  OF  THE  DATA  TO  BE  PLOTTED  IS  IN  THE 
C  FORM  OF  REAL  NUMBERS,  SO  ALL  OF  THE  VARIABLES  ARE  MADE 
C  REAL.  D  IS  USED  FOR  THE  FREQUENCY  DATA,  G  IS  USED  FOR  THE 
C  PISTON  VELOCITY  DATA,  H  IS  USED  FOR  THE  ACOUSTIC  POWER 
C  DATA,  0  IS  USED  FOR  THE  ELECTROACOUSTIC  EFFICIENCY  DATA 
C  AND,  T  IS  USED  FOR  THE  VOLTAGE  DATA.  THE  XMAX  AND  YMAX 
C  VARIABLES  ARE  INITIALIZED  TO  ZERO  AND  ARE  USED  TO  CALCU- 
C  LATE  THE  MAX  VALUE  OF  EACH  OF  THE  PARAMETERS  AS  THEY  ARE 
C  READ  INTO  THE  ARRAYS.  THIS  COMPLETES  THE  INTRODUCTION  AND 
C  EXPLAINS  THE  FIRST  FIVE  LINES  OF  THE  PROGRAM. 


DIMENSION  D{300) ,G(300) ,H(300) ,0(300) ,T{300) 

CHARACTER*80  NEW1$,DUMMY$ 

REAL  D,G,H,0,T,M 

DATA  XMAX , YMAXl , YMAX2 , YMAX3 , YMAX4/0 .,0.,0.,0.,0./ 


C  THE  NEXT  SEVEN  LINES  OF  THE  PROGRAM  RENAME  MY  DATA 
C  FILE  TO  THE  FORMAT  NECESSARY  FOR  USE  IN  DISSPLA,  OPEN  THE 
C  FILE,  READ  THE  MASS  AND  VOLUME  OUT  OF  THE  FILE  AND,  READ 
C  THE  HEADINGS  OUT  OF  THE  FILE. 


NEW1$=' RENAME  '//'  FRIDGE  OUTPUT  Al'//'  FILE  XYZl  Al 
CALL  EXCMS{NEW1$) 
OPEN  {69,FILE='XYZ1 ' ) 
READ  {69.'(A80)')  DUMMY$ 
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READ  (69,  DM,  VOL 
READ  (69,'(A80)')  DUMMY$ 
1  FORMAT  (E13.4,E13.4) 


C  THE  FOLLOWING  DO  LOOP  READS  THE  PLOTTING  DATA  INTO  THE 
C  ABOVE  SPECIFIED  ARRAYS  FOR  USE  IN  THE  DISSPLA  SUBROUTINES. 
C  AND  ALSO  CALCULATES  THE  MAXIMUM  VALUE  OF  EACH  OF  THE 
C  ARRAYS  FOR  SETTING  UP  THE  SCALES  OF  THE  PLOT  AXES.  THEN 
C  THE  DATA  FILE  IS  CLOSED.  THESE  VALUES  ARE  LEFT  AS  REAL 
C  NUMBERS  VICE  INTEGERS  BECAUSE  THE  DATA  FROM  THE  PROGRAM 
C  CAN  VARY  SIGNIFICANTLY  WITH  DIFFERENT  INPUT  VALUES,  AND  IF 
C  THE  MAX  VALUE  IS  LESS  THAN  10  AN  INTEGER  REQUIREMENT  GIVES 
C  A  PLOTTING  ERROR. 


DO  20  J=l ,300 

READ  (69, 15,END=21)D( J) .G(J),H(J),0'"^,T(J) 
15      FORMAT  (  El  3.  5,  El  3.  5,  El  3.  5,  E 13.  5,  El  3.  5) 
IF(D(J) .GT.XMAXjTHEN 
XMAX=D( J) 
ENDIF 

I F ( G ( J ) . GT . YMAX 1 ) THEN 
YMAX1=G( J) 
ENDIF 

IF ( H ( J ) . GT . YMAX2 ) THEN 
YMAX2=H( J) 
ENDIF 

IF(0(J) .GT.YMAX3)THEN 
YMAX3=0( J) 
ENDIF 

IF ( T ( J ) . GT . YMAX4 ) THEN 
YMAX4=T(J) 
ENDIF 

20  CONTINUE 

21  CLOSE (69) 


C  THE  NEXT  THREE  LINES  TAKE  CARE  OF  ANOTHER  ODDITY  OF 
C  THE  DATA,  IN  THE  CASE  OF  THE  VOLTAGE  PLOT.  DEPENDING  ON 
C  THE  INPUT  DATA  THE  VOLTAGE  WILL  LOCK  IN  AT  8  VOLTS,  AS 
C  EXPECTED,  BUT  THIS  IS  ALSO  THE  MAX  VALUE,  SO  YOU  GET  A 
C  STRAIGHT  LINE  AT  THE  TOP  OF  THE  PLOT  MARGIN.  THE  NEXT 
C  THREE  LINES  INCREASE  THE  MAX  VALUE  TO  9  VOLTS  TO  MAKE  THIS 
C  PLOTTING  CASE  MORE  CLEAR. 


IF ( YMAX4 . LE . 9 . 0 ) THEN 

YMAX4=9.0 

ENDIF 


15: 


C  THE  NEXT  EIGHT  LINES  CALCULATE  THE  INCREMENTAL  STEPS 
C  FOR  EACH  AXIS  OF  THE  VARIOUS  PLOTS.  THE  X  AXIS  IS  THE 
C  SAME  FOR  ALL,  AND  IS  LEFT  IN  INTEGER  FORMAT  FOR  CLARITY. 
C  EACH  AXIS  IS  SET  UP  FOR  TEN  INCREMENTS.  THE  DATA  FILE  IS 
C  THEN  RENAMED  TO  THE  DESIRED  FORMAT  FOR  MY  STORAGE. 


INC=INT(XMAX/10. ) 

YINC1=YMAX1/10. 

YINC2=YMAX2/10. 

YINC3=YMAX3/10. 

YINC4=YMAX4/10. 

MAX=INT(XMAX) 

NEW1$=' RENAME  '//'  FILE  XYZl  Al'//'  FRIDGE  OUTPUT  Al ' 

CALL  EXCMS(NEW1$) 


C  THIS  IS  THE  FIRST  PLOTTING  SECTION.  IT  DOES  THE 
C  INITIAL  SET  UP  OF  ALL  OF  THE  CALLS  THAT  CAN  BE  COMMON  AND 
C  PLOTS  THE  VOLTAGE  VERSUS  FREQUENCY  PLOT.  CALL  TEK618  SETS 
C  UP  THE  REQUIRED  MACHINE  FOR  THE  PLOTTING  AND  LINKS  YOU  TO 
C  THE  PLOTTER.  CALL  BLOWUP  IS  USED  TO  MAGNIFY  THE  PLOT  FOR 
C  A  BETTER  PICTURE.  CALL  YAXANG  PUTS  THE  AXIS  NUMBERS 
C  PERPINDICULAR  TO  THE  PLOT  AXIS  TO  PREVENT  RUNNING  INTO 
C  EACH  OTHER.  CALL  PAGE  SETS  UP  THE  PAGE  SIZE  IN  INCHES. 
C  CALL  AREA2D  SETS  UP  A  TWO  DIMENSIONAL  PLOT  AREA  OF  THE 
C  SPECIFIED  SIZE  IN  INCHES.  CALL  SCMPLX  SELECTS  THE  LET- 
C  TERING  TYPE  FOR  THE  HEADING.  CALL  MESSAG  AND  CALL  REALNO 
C  ALLOW  ME  TO  PLACE  A  MESSAGE  IN  THE  PLOT  AREA  AT  A  SPECIFIC 
C  POINT  SPECIFIED  IN  INCHES  IN  THE  CALL.  CALL  XNAME  AND 
C  YNAME  ALLOW  ME  TO  LABEL  THE  AXES.  CALL  HEADIN  ALLOWS  ME 
C  TO  PLACE  A  TITLE  ON  THE  PLOT.  CALL  THKFRM  AND  FRAME  ALLOW 
C  ME  TO  PUT  FRAMES  AROUND  THE  PLOT  AREA  AND  THE  WHOLE  PLOT 
C  OF  A  SET  THICKNESS.  CALL  GRAF  SETS  UP  THE  PHYSICAL  ORIGINS 
C  AND  MAX  VALUES  AND  THE  INCREMENTS  FOR  EACH  AXIS.  CALL 
C  CURVE  TELLS  THE  MACHINE  TO  PLOT  THE  DESIRED  ARRAYS  OF  DATA 
C  POINTS.  CALL  ENDPL  ENDS  THAT  SPECIFIC  PLOT  AND  GIVES  A 
C  NEW  PAGE  FOR  THE  NEXT  PLOT. 


CALL  TEK618 

CALL  BLOWUP (1.5) 

CALL  YAXANG(O) 

CALL  PAGE (11. 8. 5) 

CALL  AREA2D(8. ,6. ) 

CALL  SCMPLX 

CALL  MESSAG( 'MASS=$' , 100,0. 1,5.5) 

CALL  REALNO (M, -2, 'ABUT' , 'ABUT' ) 

CALL  MESSAG (  ' KG$ '  , 100  ,  ' ABUT '  ,  ' ABUT '  ) 

CALL  MESSAG( ' VOLUME=$' , 100,0 .1,5.) 

CALL  REALNO { VOL , -2 , ' ABUT ' , ' ABUT ' ) 

CALL  MESSAG (  'M**3$'  , 100,  'ABUT'  ,  'ABUT'  ) 
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CALL  XNAME ( ' 

CALL  YNAME ( ' 

CALL  HEADIN( 

CALL  THKFRMC 

CALL  FRAME 

CALL  GRAF (0.0, INC, MAX 

CALL  CURVE (D,T, 300,0) 

CALL  ENDPL(O) 


FREQUENCY (IN  HERTZ )$', 100 ) 
VOLThGE  (IN  VOLTS)$' , 100) 
•VOLTAGE  VS.  FREQUENCY$ ' , 100 . 2 . 
.02) 


7.5, YINC4, YMAX4) 


1) 


THIS  SECTION  PLOTS  PISTON  VELOCITY  VERSUS  FREQUENCY 


CALL  AREA2D(8. ,6. ) 

CALL  MESSAG( 'MASS=$' , 100.0. 1,5.5) 

CALL  REALNO ( M , -2 , ' ABUT ' , ' ABUT ' ) 

CALL  MESSAG( ' KG$ ' , 100, 'ABUT* , 'ABUT' ) 

CALL  MESSAG ( ' VOLUME=$ ',100,0.1,5.) 

CALL  REALNO ( VOL , -2 , ' ABUT ' , ' ABUT ' ) 

CALL  MESSAG ( •M**3$' ,100, 'ABUT' . 'ABUT' ) 

CALL  XNAME ( 'FREQUENCY (IN  HERTZ )$', 100 ) 

CALL  YNAME(  'PISTON  VELOCITYdN  METERS/SEC  )$',  100  ) 

CALL  HEADIN( 'PISTON  VELOCITY  VS.  FREQUENCY$ ' , 100 , 2 

CALL  THKFRMl .02) 

CALL  FRAME 

CALL  GRAF ( 0 . 0 , INC , MAX ,0.0, YINCl , YMAXl ) 

CALL  CURVE (D,G, 300,0) 

CALL  ENDPL(O) 


1) 


C 


THIS  SECTION  PLOTS  ACOUSTIC  POWER  VERSUS  FREQUENCY. 


CALL  AREA2D(8. ,6. ) 

CALL  MESSAG ( ' MASS=$ ' ,100,0.1,5.5) 

CALL  REALNO (M, -2, 'ABUT' . 'ABUT' ) 

CALL  MESSAG ( ' KG$ ' , 100 , ' ABUT ' , ' ABUT ' ) 

CALL  MESSAG ( ' VOLUME=$ ' , 100,0. 1 ,5. ) 

CALL  REALNO ( VOL , -2 , ' ABUT ' , ' ABUT ' ) 

CALL  MESSAG ( 'M**3$' ,100, 'ABUT' , 'ABUT' ) 

CALL  XNAME ( 'FREQUENCY (IN  HERTZ )$', 100 ) 

CALL  YNAME ( 'ACOUSTIC  POWER  (IN  N-M/S ) $ ' , 100 ) 

CALL  HEADIN( 'ACOUSTIC  POWER  VS.  FREQUENCYS ' , 100 , 2 

CALL  THKFRMl .02) 

CALL  FRAME 

CALL  GRAF(0 .0 , INC, MAX, 0 .0 , YINC2 , YMAX2 ) 

CALL  CURVE (D,H, 300, 0) 

CALL  ENDPL(O) 


1) 


THIS  SECTION  PLOTS  POWER  RATIO  VERSUS  FREQUENCY 
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CALL  AREA2D(8. ,6. ) 

CALL  MESSAG( 'MASS=$' , 100 ,0.1,5.5) 

CALL  REALNO { M , -2 , ' ABUT ' , ' ABUT ' ) 

CALL  MESSAGC ' KG$ ' , 100, 'ABUT' , 'ABUT' ) 

CALL  MESSAG( •VOLUME=$' . 100,0. 1,5. ) 

CALL  REALNO ( VOL , -2 , ' ABUT ' , ' ABUT ' ) 

CALL  MESSAG( 'M**3$' , 100, 'ABUT' , 'ABUT' ) 

CALL  XNAME{ 'FREQUENCY (IN  HERTZ )$', 100 ) 

CALL  YNAME( 'ACOUSTIC  PWR/  ELECTRIC  PWR$',100) 

CALL  HEADIN( 'POWER  RATIO  VS.  FREQUENCYS ' , 100 , 2 . 

CALL  THKFRM( .02) 

CALL  FRAME 

CALL  GRAF (0.0, INC , MAX ,0.0, YINC3 , YMAX3 ) 

CALL  CURVE(D, 0,300,0) 

CALL  ENDPL(O) 


1) 


C      CALL   DONEPL  ENDS  THE  PLOTTING  SEQUENCE  AND  SENDS 
C  BACK  TO  CMS. 


YOU 


CALL  DONEPL 

STOP 

END 
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APPENDIX  C.   STAR  DRIVER  DIMENSIONAL  DRAWINGS 


\-^- 


ZP 


Figure  C-l . 


Assembly  Drawing  of  NIB  Driver  with 
all  Housing  Parts  and  Test  Lid. 


156 


-M^^^ 


DIA 


DIA 


hlS§DiA 


A^ 


DRIVER  HOUSING  -  RESONATOR  VIEW 


[DlmeneloriB     in     inches  ) 


Figure  C-2.   Driver  Housing--Resonator  View 
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Figure   C-3. 


Driver  Housing--Side  View  Section  AA 
Showing  Electrical  Leads,  Bolt  Holes, 
Microphone  Penetration,  and  FET 
Follower  Housing. 
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Figure  C-4 


Feed-Through  for  Electrical  Leads 
Shown  in  Figure  V-5. 
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Figure   C-5 


Driver  Housing--Side  View  Section  BB 
Showing  Pressure  Gage  and  Helium 
Fill  Ports. 
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Figure  C-6. 


DC  Pressure  Gage  Feed- 
Through  Plug  Details. 
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DRIVER  HOUSING  -  LID  VIEV 


Figure  C-7. 


Driver  Housing--Lid  View  Showing 
Capillary,  Microphone,  and  Follower, 
and  Bolt  Circles  for  Resonator 
Attachment,  Pusher  Plate,  Pressure 
Lid  and  GAS  Can  Lid. 
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PUSHER    PLATE    (LID    VIEW) 
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Figure  C-8.      Aluminum   Pusher   Plate 
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PRESSURE    LID    (LID    VIEV) 


( Dimensions  in  inches  ) 


Figure  C-9. 


Pressure  Lid  with  Lead 
0-ring  Groove. 
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Figure  C-10.   Pressure  and  Helium  Leak  Test  Lid 
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Figure   C-11.       Reducer   Cone 
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APPENDIX  D.   STAR  DRIVER  LIST  OF  MATERIALS 

The  driver  consists  of  a  modified  JBL  model  2450J  neo- 
dymium-iron-boron  magnet  with  iron  pole  pieces  on  either 
side.  The  diaphragm  is  titanium  and  the  voice  coil  is 
aluminum . 

The  housing.  pusher  plate,  pressure  lid,  and  reducer 
cone  will  all  be  machined  from  6061-T6  aluminum.  All  open- 
ings in  the  driver  housing  will  be  sealed  with  1/16"  or 
1/32"  diameter  lead  o-rings.  The  electrical  feed-throughs 
will  be  sealed  with  Emerson  &  Cuming's  STYCAST  <"r'*'>  2350FT 
epoxy.  All  wire  leads  are  copper.  The  steel  bolts  connecting 
the  resonator  to  the  driver  are  size  1/4-20.  All  other  steel 
bolts  are  size  10-32. 
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